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FINAL  REPORT  ON  A  CALCULATIONAL  PARAMETER 
STUDY  OF  SOILS  TYPICAL  OF  SOME  ESSEX  I  CRATERING  SITES 


Abstract 


The  one-dimensional  computer  calcula¬ 
tions  described  in  this  report  were  per¬ 
formed  to  simulate  stress-wave  propagation 
and  kinetic  energy  transfer  associated  with 
subsurface  cratering  detonations  in  soils. 

A  hypothetical  20-ton-yield  nuclear  explo¬ 
sive  was  assumed  as  the  energy  source, 
surrounded  by  a  single  soil  material. 
Various  soil  descriptions  were  selected  in 


order  to  systematically  study  the  range  of 
soil  response  to  the  nuclear  detonation. 
The  soils  were  representative  of  the  lay¬ 
ered  mixtures  of  sand  and  clay  found  at 
the  ESSEX  high-explosive  cratering  sites 
near  Ft.  Polk,  Louisiana.  Soil  properties 
analyzed  in  this  study  include  water  sat¬ 
uration,  bulk  density,  failure  envelope, 
and  low-pressure  bulk  modulus. 


Introduction 


The  Earth  Sciences  (K)  Division  of  the 

Lawrence  Livermore  Laboratory  (LLL)  has 

conducted  a  calculational  parameter  study 

* 

for  the  Explosive  Excavation  Division  of 
the  Weapons  Effects  Laboratory  of  the  U.S. 
Army  Engineer  Waterways  Experiment  Station. 
The  study  began  a  .systematic  investigation 
into  the  influence  of  soil  properties  on 
nuclear  cratering  and  ground  motion. 

Work  performed  during  the  study  included 
investigating  the  static  properties  of 
soils,  determining  the  relevant  range  of 
values  of  those  properties  to  be  used  in 


the  study,  developing  a  constitutive  rela¬ 
tions  model  that  could  derive  representa¬ 
tive  dynamic  properties  from  given  static 
properties  for  hypothetical  soils  spanning 
the  ranges  of  interest  in  the  study,  and 
calculating  the  response  of  each  hypothet¬ 
ical  soil  to  a  nuclear  detonation  in  one 
dimension. 

This  report  reviews  the  work  performed, 
with  emphasis  on  the  results  of  the  calcu¬ 
lations.  For  some  topics,  more  detailed 
information  can  be  found  in  the  two  prog- 
ress  reports. 


Objective  and  Concept 


OBJECTIVE 

The  calculational  parameter  study  de¬ 
scribed  here  is  part  of  a  long-range 


Formerly  known  as  the  Explosive  Excavation 
Research  Laboratory  (EERL)  of  the  U.S.  Army 
Engineer  Waterways  Experiment  Station  (WES) ; 
located  at  LLL,  Livermore,  California. 


parameter  study  proposed  by  the  Explosive 
3 

Excavation  Division.  The  overall  objec¬ 
tive  of  this  study  is  to  examine  the 
phenomenology  of  nuclear  cratering  and  to 
present  the  results  in  a  form  directly 
usable  in  an  improved  U.S.  Army  field 
handbook  predicting  various  effects  of 
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Fig.  1.  Concept  of  total  nuclear  effects  study.  The  study  will  produce  information  for  a 
handbook  that  will  predict  nuclear  phenomenology  given  the  geology,  yield,  and 
emplacement. 


nuclear  weapons.  Generally  these  effects 
are  influenced  by  the  site  geology,  the 
explosive  yield,  and  the  emplacement 
(i.e.,  both  the  configuration  and  the 
depth  of  burial).  The  effects  of  greatest 
interest  are  ground  motion,  crater  size, 
and  ancillary  cratering-related  effects. 

The  concept  of  the  total  study  is  shown 
schematically  in  Fig.  1. 

A  complete  investigation  into  the 
effects  of  site  geology  alone  requires  more 
work  than  is  entailed  in  this  year’s  calcu-? 
lational  parameter  study.  Thus  it  is  only 
a  beginning  step  toward  accomplishment  of 
the  long-range  objectives  of  the  study  and 
is  not  expected  to  answer  all  questions 
about  cratering  effects.  The  primary  goal 
of  this  year’s  study  was  to  identify  the 
material  properties  that  have  greatest 
influence  on  energy  coupling  and  stress- 
wave  propagation  resulting  from  a  buried 
nuclear  detonation. 


CONCEPT 

The  concept  of  the  calculational  parame¬ 
ter  study  is  shown  schematically  in  Fig.  2. 
Material  properties  were  separated  into  two 
groups:  bulk  (or  initial)  properties  and 
constitutive  relations.  Bulk  properties 
relate  to  the  undisturbed  soil  only,  where¬ 
as  constitutive  relations  describe  the 
effects  of  pressure  on  the  material.  A  com¬ 
prehensive  (although  not  exhaustive)  list 
of  bulk  properties  was  compilea.  Kelation- 
ships  between  these  properties  were  derived, 
and  a  set  of  independent  variables  was 
chosen. 

Since  this  is  a  follow-on  study  designed 
to  clarify  and  expand  the  results  of  the 
eight  ESSEX  cratering  tests,  conducted 
in  the  Peason  Ridge  area  of  Fort  Polk, 


*ESSEX  is  an  acronym  formed  from  (effects 
of  ^ub^urface  explosions . 
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Louisiana,  *  the  study  was  limited  to 
soils  characteristic  of  that  site.  The 
Soils  and  Pavements  Laboratory  (S  &  PL)  of 
WES  conducted,  for  ESSEX,  an  extensive 
investigation  of  the  soils  at  that  test 
site.  The  investigation  produced  13  ideal¬ 
ized  units  with  material  properties  typical 
of  subsurface  materials  encountered  at  the 
eight  locations  studied.^  These  units  were 
examined  and  used  to  infer  appropriate 
limits  for  the  study’s  independent 
variables. 

Since  a  knowledge  of  the  constitutive 
relations  was  needed  in  order  to  calculate 
the  effect  of  a  nuclear  detonation  in  a 
given  material,  a  model  was  developed  to 
supply  representative  constitutive 
relations  given  the  values  of  the 
independent  bulk  properties.  The  model’s 
applicability  was  established  by  comparing 
model-generated  constitutive  relations  for 
the  13  Fort  Polk  units  with  those  recom¬ 
mended  by  S  &  PL  (based  on  laboratory 
tests).  The  model  was  then  used  to  develop 
constitutive  relations  for  hypothetical 
soils  with  bulk  properties  spanning  the 
ranges  proposed  for  the  study. 


One-dimensional,  spherical  code  calcula¬ 
tions  were  performed  for  each  hypothetical 
soil  by  means  of  the  SOC74  code,  which  is  a 
time-dependent,  Lagrangian,  finite- 
difference  computer  code^  that  can 
simulate  stress-wave  propagation  through 
solid  materials  with  generalized 
characteristics. 

Stress  wave  propagation  away  from  the 
nuclear  source  and  energy  coupling  into 
the  soil  were  calculated  as  a  function  of 
time.  The  results  were  plotted,  analyzed, 
and  compared  with  other  calculational 
results.  Sensitivity  of  the  results  to 
individual  bulk  properties  could  then  be 
investigated  systematically. 

The  general  method  became  one  of 
performing  calculations  on  hypothetical 
materials,  and  varying  the  independent 
bulk  parameters  within  the  set  limits. 

When  the  results  proved  sensitive  to  a 
particular  parameter,  more  calculations 
were  performed  using  additional  values  of 
that  bulk  property  within  the  limiting 
range.  This  process  was  continued  until 
adequate  curves  could  be  drawn  showing  the 
sensitivity  of  the  calculational  results  to 
that  bulk  property. 


Material  Properties  and  Constitutive  Model 


GENERAL  CONSIDERATIONS 

In  this  study  material  properties  have 
been  considered  under  two  categories:  bulk 
properties  and  constitutive  relations. 

Bulk  properties  are  static  characteristics 
of  the  material  in  its  undisturbed  state 
(initial  condition)  only,  whereas  constitu¬ 
tive  relations  define  the  effects  of 
pressure  on  the  materials. 


BULK  MATERIAL  PROPERTIES 

A  review  of  bulk  parameters,  which 
quantitatively  relate  information  about 
bulk  properties,  was  performed  for  the 
purpose  of  defining  independent  parameters. 
A  comprehensive  but  not  exhaustive  list 
of  bulk  parameters  was  compiled  and  is 
shown  in  Table  1.  Mathematical  relation¬ 
ships  between  these  parameters  were  derived 
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Table  1.  Definition  of  bulk  material  parameters. 


Parameters 

S3nnbol  used 
in  study 

It 

CO 

Description 

Densities 

PQ 

Composition  parameters^ 

Bulk  density  (sample  mass/sample  volume) 

Ps 

Dry  density  (mass  of  solids/sample  volume) 

Pw 

Mass  of  water  per  unit  total  volume  (water  mass/ 
sample  volume) 

Pg 

G 

s 

Grain  density  (mass  of  solids/volume  of  solids) 

Volumes 

e 

e 

Void  ratio  (volume  of  voids/volume  of  solids) 

e 

w 

— 

Water  ratio  (volume  of  water/volume  of  solids) 

-*>0 

-  — 

Total  porosity  (volume  of  voids/sample  volume) 

V 

a 

Air-filled  porosity  (air  volume /sample  volume) 

V 

w 

Water-filled  porosity  (water  volume/ sample  volume) 

‘f’s 

V 

s 

Solid  volume  (volume  of  solids/sample  volume) 

Saturations 

S 

Water  saturation  (water  volume/volume  of  voids) 

s 

a 

— 

Air  saturation  (air  volume/volume  of  voids) 

Weights 

W 

w 

Water  content  (weight  of  water /dry  weight  of  sample) 

Z 

— 

Water  content  (weight  of  water /wet  weight  of  sample) 

Elastic  parameters^ 

K 

Bulk  modulus 

G 

Shear  modulus 

*V 

Poisson’s  ratio 

A 

Lame  elastic  parameter 

E 

Young’s  modulus 

B 

Constrained  modulus 

^0 

Stress  ratio  in  uniaxial  strain 

Three  independent  parameters. 
^Two  independent  parameters. 


and  are  given  in  Table  2.  These  parameters 
fall  naturally  into  two  groups:  composition 
parameters  and  elastic  parameters.  Exami¬ 
nation  of  the  relationships  between  each 
group’s  members  indicates  a  total  of  five 


independent  bulk  parameters.  Of  these 
five  independent  parameters,  three  must 
be  composition  and  two  elastic.  Composi¬ 
tion  parameters  chosen  as  independent  must 
include  one  density  and  one  volume  ratio, 
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Table  2,  Mathematical  relationships  among  parameters. 


Composition  parameters  (using  p^,  p^,  and  as  a  base) 
5^  (Note:  the  density 

of  water  is  unity) 


Pq  =  Pg(l  -  (Jiq)  +  (Note:  the  density  <j)^  =  (1)q(1  -  S^)  (Note:  <(>3  =  1 


P,  -  PgCl  -  ♦»> 


'^O^W 


‘■g-Og 


e  -  *^/a  -  *a>  ■ 


6  =  (h^S 


=  1  - 


S  =  1  -  S„  (Note:  S,,  +  S  =  1) 
a  W  W  a 


%  ■  ■  ♦oV«  -  ■  <'»o  -  -  *0*  “  ■  '»oV>''g<i  -  *0)1  ■ 

(P„  "  Pn)/(P„  "  S  ) 


g  0  g  w 


Z  -  ♦oVff's'"  ■  ♦o>  *  Vw'  -  »yPo 


K  =  K 


Elastic  parameters  (using  K  and  V  as  a  base) 

E  =  3K(1  -  2v) 


G  =  3K(1  -  2v)/2(l  +  V) 


B  =  3K(1  -  V)/(l  +  V) 


V  =  V 

X  =  3Kv/a  +  V) 


or  saturation.  The  remaining  composition 
parameter  and  the  two  independent  elastic 
parameters  can  then  be  selected  from  their 
respective  groups,  at  the  investigator’s 
convenience. 

Since  the  study  was  limited  to  materials 
typical  of  Fort  Polk,  the  value  of  each 
bulk  parameter  was  calculated  for  all  13 
Fort  Polk  units.  These  are  shown  in 
Table  3;  the  ranges  for  each  of  these 
parameters  will  be  discussed  in  the  next 
section.  Values  of  the  elastic  parameters, 
other  than  the  bulk  modulus  K  and  the 
limiting  value  of  Poisson’s  ratio  V,  are 
AOt  shown. 


k^  =  V/(1  -  V) 


CHOOSING  INDEPENDENT  BULK  PARAMETERS 

The  three  composition  parameters  cuoocii. 

as  independent  during  the  study  were  grain 

density  p^,  bulk  density  p^,  and  water 

saturation  S  .  These  were  chosen  primarily 
w 

for  practical  reasons.  Grain  density  can 
be  assumed  to  be  constant  (p  =2.67  Mg/m  ), 
since  it  shows  only  a  slight  variation  at 
Fort  Polk  (2.66  to  2.70  Mg/m^,  see  Table  3). 
It  was  felt  that  bulk  density  would  be 
more  easily  measured  under  field  condi¬ 
tions  than  most  other  pa'rameters.  Water 
saturation  was  used  because  previous 
experience  on  rocks  indicates  that  this 
factor  has  a  strong  influence  on  both 
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ground  motion  and  cratering  in  the  range 
of  interest  (0.9  to  1.0).^^’^^ 

Bulk  modulus  K  and  Poisson’s  ratio  V 
were  chosen  as  the  independent  elastic 
parameters.  Bulk  modulus  was  chosen 
because  it  can  be  estimated  from  a 
measurement  of  the  speed  of  sound  if  the 
initial  shear  modulus  is  small  enough  to 
be  neglected.  Poisson’s  ratio  increases 
very  rapidly  with  confining  pressure  at 
very  low  pressures  (-.01  GPa)  from  its 
initial  value,  shown  for  each  unit  in 
Table  3,  to  a  value  of  0.46  to  0.48 
(very  close  to  the  limiting  value  at  high 
pressures  of  0.5)  for  almost  all  Fort  Polk 
units.  Hence,  instead  of  using  the  range 
of  initial  values  of  V  from  the  Fort  Polk 
data,  a  constant  value,  0.48,  was  used. 

Therefore  two  of  the  five  independent 
parameters  were  held  constant,  while  three 
were  varied.  Limits  of  variation  for  these 
three  variables  were  chosen  so  as  to 
completely  cover  their  ranges  at  Fort  Polk. 
Table  4  lists  the  values  of  the  two  inde¬ 
pendent  bulk  parameters  held  constant  and 
the  ranges  of  values  of  the  remaining  three. 

By  using  the  values  of  the  independent 
composition  parameters  (p^,  p^,  and  S^) 

Table  4.  Range  of  values  of  the  five  inde¬ 
pendent  bulk  parameters  consid¬ 
ered  in  the  study 

Independent  bulk  parameters  held  constant 

Grain  density  p^:  2,67  Mg/m 

Poisson’s  ratio  v:  0.48 

Independent  bulk  parameters  varied 

Bulk  density  p^:  1.7  to  2,1  Mg/m 

Water  saturation  S  :  0,8  to  1.0 

w 

Initial  bulk  modulus  K:  2.5  to  7.5  GPa 


shown  in  Table  4,  the  ranges  of  each 
dependent  composition  parameter  exhibited 
by  the  13  Fort  Polk  units  are  spanned 
more  than  adequately.  This  is  shown  in 
Table  5,  which  compares  the  ranges  of  the 
dependent  composition  parameters  calculated 
from  Table  4  with  the  ranges  obtained 
directly  from  Table  3.  Not  only  the  limits 
but  also  the  distribution  of  the  Fort  Polk 
data  within  those  limits  is  important  for 
the  study.  Figure  3  shows  this  distribu¬ 
tion  by  plotting  p^  versus  for  the  13 
Fort  Polk  units.  The  data  are  randomly 
scattered  within  an  envelope  approximately 

Table  5.  Ranges  for  the  dependent 
composition  parameters. 


Composition 

parameter 

Range 
Calculated 
from  Table  4 

Range  from  13 
Fort  Polk  units 
(Table  3) 

Ps 

1.12-1.86 

1.2417-1.6546 

Pw 

(Mg/m^) 

0.24-0.58 

0.3230-0.5318 

e 

— 

0.6137-1.1744 

<^0 

0.30-0.58 

0.3803-0.5401 

0.0-0.12 

0.0-0.069 

K 

0.24-0.58 

0.3230-0.5318 

e 

w 

0.5317-1.136 

0.42-0.70 

0.4599-0.6197 

S 

a 

0.0-0.20 

0.0-0.17 

W 

~ 

0.1999-0.4207 

z 

0.11-0.34 

0.167-0.296 

-8- 


limited  by  a  range  from  1.75  to  2.05  and 

an  S  range  of  0.9  to  1.0,  excluding  only 
w 

unit  1.  The  parameter  study  covers  a  some¬ 
what  larger  area,  with  ranging  from  1.7 

to  2.1  and  S  ranging  from  0.8  to  1.0. 
w 

This  was  done  to  ensure  complete  coverage 
of  the  dependent  variables  over  their 
respective  ranges  at  Fort  Polk  and  account 
for  the  possibility  that  the  composition 
of  all  materials  at  the  site  might  exhibit 
a  wider  variation  than  that  determined 
from  the  13  idealized  units. 

For  the  elastic  parameters,  the  range 
of  values  of  initial  bulk  modulus  (2.5  to 
7.5  GPa)  does  not  span  the  entire  range  of 
values  found  in  the  Fort  Polk  data  (0.172 
to  7.31  GPa),  as  seen  from  Table  3. 
Experience  has  shown  that  the  effect  of 
the  elastic  parameters  on  nuclear  crater¬ 
ing  and  close-in  ground  motion  (3  to  15  m 
for  an  84-GJ  nuclear  detonation)  was 
minuscule.  A  few  calculations  were  per¬ 
formed  at  different  values  of  K  to  prove 
this  point,  and  no  further  analysis  of 
the  remainder  of  the  elastic  parameters 
shown  in  Table  1  was  performed. 

CONSTITUTIVE  MODEL  FOR  SOILS 

A  constitutive  model  was  developed  to 
provide  constitutive  relations  for  hypo¬ 
thetical  soils  having  bulk  properties  in 
the  ranges  discussed  in  the  preceding  sec¬ 
tion.  This  model  consists  of  two  parts: 
a  computer  model  that  predicts  compressi¬ 
bility  and  an  empirical  model  that  gives 
failure  information.  This  section  dis¬ 
cusses  briefly  each  of  these  models.  More 

detailed  information  is  presented  in  a 

,  2 

previous  progress  report. 

Since  the  S0C7A  computer  code  was  used 
for  the  calculations,  the  model  was  con- 


0.7  0.8  0.9  1.0 


Water  saturation 

Fig.  3.  Distribution  of  Fort  Polk  units  in 
space.  The  number  beside 
each  data  point  is  the  Fort  Polk 
unit  number. 

structed  so  as  to  output  constitutive  re¬ 
lations  compatible  with  that  code.  The 
SOC74  code  provides  for  a  great  deal  of 
flexibility  in  the  modeling  of  material 
behavior.  Regimes  of  response  that  can 
be  modeled  include  linear  elastic  or 
incrementally  elastic  compression  hyster- 
etic  compaction,  brittle  failure,  ductile 
flow  (elastic-plastic  failure)  tensile 
failure,  liquefaction  or  vaporization  of  a 
solid,  vaporization  of  only  the  water 
component  in  a  solid,  and  gaseous  be¬ 
havior.  Techniques  are  available  to  sim¬ 
ulate  certain  types  of  rate-dependent 
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Pressure  P 


behavior.  Since  insufficient  data  are 
available  to  characterize  rate-dependent 
effects  in  the  soils  being  considered,  the 
constitutive  relations  for  all  hypothetical 
soils  used  in  the  study  were  assumed  to  be 
rate- independent.  The  models  used  in  the 
SOC74  code  also  assume  that  materials  are 
isotropic  and  homogeneous. 

The  SOC74  code  uses  tabular  stress- 
strain  relationships  to  specify  the  re¬ 
sponse  of  solid  media  under  hydrostatic 
compression.  The  user  enters  the  tables, 
which  give  the  pressure  P  as  a  function  of 
the  excess  compression  y  (y  =  p/Pq 
where  is  the  initial  density  of  the 
material  and  p  is  its  density  at  the  pres¬ 
sure  P) .  The  general  form  of  the  relation 
expected  by  the  code  is  shown  in  Fig.  4. 

Two  curves  are  required  for  each  material: 
a  "virgin  loading"  curve  and  a  "completely 


Fig.  4.  General  form  of  the  compressibility 
curve  for  the  SOC74  code. 


t 


Fig.  5.  General  form  of  the  failure  curve 
for  the  SOC74  code. 


crushed"  unloading  (and  reloading)  curve. 
The  virgin  loading  curve  is  assumed  to  be 
reversibly  elastic  until  the  transition 
pressure  P^  is  exceeded;  once  exceeded, 
irreversible  compaction  (hysteretic  be¬ 
havior)  is  allowed.  The  loading  and  un¬ 
loading  curves  are  assumed  to  merge  into 
a  single  curve  at  and  above  a  specified 
merge  pressure  P^.  All  air  void  space  is 
assumed  to  be  irreversibly  removed  on 

loading  between  P^  and  P  .  Partial  com- 
t  m 

paction  is  allowed  on  loading  to  a  maximum 

pressure  between  P  and  P  .  No  further 
t  m 

hysteretic  compaction  is  allowed  at 

pressures  above  P  . 

m 

The  failure  criterion  in  SOC74  is  repre¬ 
sented  by  a  table  giving  permissible  shear 
stress  T  versus  confining  pressure  P. 

This  specification  limits  the  deviator ic 
stress  (shear  stress)  that  a  material  can 
support.  Brittle  failure  and  strength 
reduction  may  be  simulated  by  specifying 
two  strength  curves  for  a  material:  a 
"virgin"  curve  and  a  "completely  failed" 
curve.  The  material  strength  is  gradually 
reduced  from  the  virgin  curve  to  the 
completely  failed  curve  as  damage  to  the 
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material  increases.  The  general  form  of 
the  relation  expected  by  SOC74  is  shown 
in  Fig.  5. 

Compressibility  Model 

The  soil-compressibility  model  is  an 

extension  of  the  model  developed  by 
13 

Butkovich  for  certain  rock  types.  Input 
numbers  for  the  computer  model  are  shown 
in  Table  6.  Three  of  the  first  four  input 
numbers,  p^,  S^,  and  K,  will  be  recognized 
as  the  independent  bulk  parameters  chosen 
for  variation  in  the  study  (see  Table  4). 

Values  for  P_  and  P  are  dependent 
1  m 

somewhat  on  the  failure  model  used;  input 
for  these  quantities  is  discussed  in  the 
section  that  follows.  Finally,  the  weight 
fractions  of  clay  and  sand  must  be  pro¬ 
vided  so  that  proper  weighting  of  the 
loading  Hugoniots  of  sand  and  clay  can  be 
computed,  as  discussed  in  Ref.  2.  These 
weight  fractions  were  held  constant  (ZCLA  = 
0.25  and  ZSIL  =  0.75)  in  the  study. 

The  model  appears  to  be  quite  adequate 

for  predicting  compressibility  for  soils 

typical  of  Fort  Polk.  In  a  recent  prog- 
2 

gress  report,  compressibility  curves 
generated  using  the  bulk  properties  of  the 

Table  6.  Input  required  for  soil-compres¬ 
sibility  model 

Symbol 

3 

Initial  bulk  density  (Mg/m  ) 

^  3 

pg  Grain  density  (2.67  Mg/m  ) 

S  Water  saturation 

w 

K  Initial  bulk  modulus  [GPa  x  100  (Mbar) 
Transition  pressure  [GPa  x  IQO  (Mbar)] 

P^  Merge  pressure  [GPa  x  100  (Mbar)] 

ZCLA  Weight  fraction  of  clay  (0.25) 

ZSIL  Weight  fraction  of  sand  (0.75) _ 


Fig.  6.  Comparison  of  the  model-generated 
compressibility  curve  using  the 
bulk  properties  of  Fort  Polk  unit  5 
soil  with  the  WES-recommended  curve. 


13  Fort  Polk  units  (see  Table  3)  were  com¬ 
pared  to  the  WES-recommended  curves  (S  & 
PL),  which  were  based  on  laboratory  tests. 
For  almost  all  units,  adequate  agreement 
was  observed.  Figure  6  shows  this  com¬ 
parison  at  low  pressures  (0  to  0.04  GPa) 
for  unit  5.  This  unit  layer  was  recom¬ 
mended  at  or  near  shot  depth  on  several  of 
the  ESSEX  I  high-explosive  experiments.^ 
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Failure  Model  and  Strength  Parameters 
As  discussed  earlier  in  this  section, 
SOC74  requires  a  failure  curve  of  the  form 
shown  in  Fig.  5.  Moreover,  the  compressi¬ 
bility  model  requires  the  input  of  two 
strength  parameters,  and  P^,  So  far 
three  independent  bulk  parameters  have  been 
chosen  for  variation.  If  in  addition  the 
failure  curve  and  P..  and  P^  were  also 
independently  varied,  the  number  of  in¬ 
dependent  variables  in  the  study  would 
become  prohibitively  large.  This  section 
discusses  how  the  failure  curves  for  the 

hypothetical  soils  were  chosen  and  how  P,, 

M 

and  P^  were  tied  to  those  curves,  thereby 
reducing  the  number  of  independent  varia¬ 
bles  from  six  to  four. 

2 

In  the  progress  report  an  attempt  was 
made  to  correlate  the  measured  maximum 
shear  strengths  of  the  Fort  Polk  soils 
with  their  bulk  properties,  in  particular, 
water  content.  Such  a  correlation  was 


found  by  Butkovich  to  exist  for  some 
rock  types.  However,  the  results  found 
for  Fort  Polk  soils  were  inconclusive. 

In  order  to  establish  approximate  fail¬ 
ure  criteria  for  the  hypothetical  soils, 
then,  an  empirical  model  was  developed. 

This  was  based  on  the  WES-recommended  fail¬ 
ure  curves  for  the  13  Fort  Polk  units^  and 

14 

for  an  additional  unit,  unit  14.  This 
latter  was  identified  by  S  &  PL  as  a  "weak 
rock,"  having  strength  characteristics 
exceeding  those  of  the  other  13  units. 

These  units  were  divided  into  four  groups; 
uncement ed-satur a ted  (units  6,  11,  and  12), 
uncement ed-un saturated  (units  1,  2,  4,  5, 

7,  8,  and  10),  cement ed-unsaturated  (units 
9  and  13),  and  weak  rock  (unit  14).  Maxi¬ 
mum  and  minimum  typical  (or  fitted)  curves 
were  then  drawn  for  each  of  the  four 
groups.  These  sets  of  typical  curves  were 
then  used  as  another  independent  variable. 


Fig.  7.  Failure  curves  for  uncemented- saturated  soils  (SM  and  SX) , 
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Fig.  8.  Failure  curves  for  uncemented-unsaturated  soils  (UM  and  UX) . 


Typical  curves  for  the  uncemented- 
saturated  group  (SM  and  SX)  are  shown  in 
Fig,  7.  In  this  case  they  correspond  to 
the  actual  maximum  and  minimum  curves 
recommended  by  WES.  Fig.  8  shows  typical 
curves  for  the  uncemented-unsaturated 
group  (UM  and  UX) .  They  also  correspond 
closely  with  the  maximum  and  minimum 
curves  recommended  by  WES.  The  unsaturated- 
cemented  (CM  and  CX)  curves,  as  Fig.  9 
shows,  do  not  follow  the  data  as  closely. 

In  particular,  the  plateau  or  step  that 
appears  in  the  data  for  unit  9  was  not 
included  in  the  CX  curve.  This  was  because 
this  region  is  not  reached  on  initial  load¬ 
ing.  (Use  of  a  constant  Poisson's  ratio 
does  not  allow  failure  to  occur  below  . 
Since  the  step  behavior  is  thought  to  be 
associated  with  pore  collapse,  it  should 
not  be  present  during  unloading  after  load¬ 
ing  to  higher  pressures.  Thus  the  same 


curve,  CX,  was  used  for  both  loading  and 
unloading.  Finally,  Figure  10  shows  the 
typical  curves  for  weak  rock  (RM  and  RX) . 

The  RX  curve  follows  the  unit  14  curve 
closely.  The  RM  curve  was  chosen  so  as  to 
have  a  maximum  shear  strength  exactly 
half  that  of  the  RX  curve,  10  MPa. 

For  the  study,  the  merge  pressure  P^^  was 
chosen  to  be  a  function  of  the  maximum 
shear  strength  Figure  11  shows  P^^ 

versus  t  for  the  unsaturated  Fort  Polk 

X 

units.  The  values  show  considerable 
scatter.  Hence,  somewhat  arbitrarily,  P^ 
(in  meg'apascals)  was  related  linearly  to  T^: 

P^  =  5t  +  9.5. 

M  X 

This  line  falls  in  the  middle  of  the  data 
shown  in  Fig.  11  and  hence  is  representa¬ 
tive  of  Fort  Polk  soils. 

The  transition  pressures  P^  were  the 
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Shear  stress  t  (MPa) 


3.0 


Fig.  9.  Failure  curves  for  unsaturated-cemented  soils  (CM  and  CX) . 


Fig.  10,  Failure  curves  for  weak  rock  (RM  and  RX) 


Table  7.  Pt.  and  used  for  hypo¬ 

thetical  soils  of  the  Fort  Polk 
type. 


X 

Prp 

Failure 

X 

M 

T 

Model 

Curve 

(MPa) 

(MPa) 

(MPa) 

Uncemented, 

SM 

0.17 

0.31 

0.01 

saturated 

SX 

0.23 

0.50 

0.01 

Uncemented, 

UM 

0.17 

10.4 

0.1 

unsaturated 

UX 

1,03 

14.7 

0.1 

Cemented, 

CM 

2,40 

21.5 

5.0 

unsaturated 

cx 

3.10 

25.1 

5.0 

Weak  rock 

RM 

5.0 

34.5 

8.0 

RX 

10.0 

60.0 

8.0 

same  for  each  set  of  curves  for  all  four 
failure  models.  Their  values  were  derived 
by  averaging  the  values  given  for  each 
group  of  Fort  Polk  units. 

Table  7  summarizes  P^^,  and  P^  for 
each  of  the  eight  failure  curves. 


Maximum  shear  strength  (MPa) 


Fig.  11.  Merge  pressure  Pj^  vs  maximum  shear 
strength  for  Fort  Polk  units. 


Calculations  Performed 


ZONING,  MATERIAL  BOUNDARIES  AND  INITIAL 
CONDITIONS 

A  hypothetical  nuclear  source  with  an 
energy  yield  of  84  GJ  (0.02  kt)  was  used 
for  all  calculations.  It  was  represented 
by  a  sphere  of  iron  gas  initially  having  a 
0.27-m  radius,  a  370-GPa  pressure,  and 

o 

a  1.5-Mg/m  density.  The  equation  of  state 
of  the  iron  gas  used  in  SOC74  was  developed 
by  Chapin  and  Butkovich.^^ 

The  nuclear  source  was  surrounded  by  a 
spherical  shell  of  hypothetical  soil  extend¬ 
ing  from  the  source  radius  to  a  radius  of 


200  m,  which  is  effectively  infinite  for 
these  problems.  Thus  there  were  no  reflec¬ 
tions  from  the  outer  boundary  during  the 
20-msec  simulation. 

Both  source  and  soil  were  divided  into 
zones  consisting  of  concentric  spherical 
shells.  The  source  contained  20  zones  of 
equal  thickness,  while  the  soil  contained 
505  zones  varying  in  thickness  according  to 
a  geometrical  progression  from  0.0136  m  at 
the  edge  of  the  source  to  about  2  m  at  the 
outer  edge  of  the  problem  (200  m) .  The 
initial  state  of  the  problems,  including 
zoning  and  initial  source  conditions,  is 
shown  diagrammatically  in  Fig.  12. 
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Material 

boundaries 


Zoning 


BULK  PROPERTY  AND  CONSTITUTIVE-RELATION 
INPUT  FOR  INDIVIDUAL  CALCULATIONS 


Fig.  12.  Diagram  of  zoning  and  material 
boundaries  for  SOC74  hypotheti¬ 
cal  nuclear  problems. 


Results  from  25  SOC74  calculations  are 
included  in  this  report.  For  convenience 
in  referring  to  individual  calculations, 
mnemonics  were  developed.  The  mnemonic  is 
of  the  following  form: 


where  F  gives  the  letter  designation  of 
the  failure  model  used  (e.g.,  UM  or  CX) . 
The  middle  part,  p^/S^,  gives  the  initial 
values  of  (Mg/m^)  and  S^  used.  There 
were  only  two  calculations  with  values  of 
K  that  differed  from  5.0  GPa.  For  these 
two  calculations  the  value  of  K  is  in¬ 
cluded  in  parentheses  to  the  right  of  the 
mnemonic.  (See  Table  8,  p.  l8.)' 

The  compressibility  curves  on  loading 
and  unloading,  for  each  of  the  25  calcula¬ 
tions,  are  included  in  Appendix  A.  This 


Failure  Model  [maximum  shear  stress  (MPa)] 
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Fig.  13.  Summary  of  calculations  performed  in  the  parameter  study. 
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inforaation,  together  with  the  appropriate 
failure  curve  (Figs.  7,  8,  9,  or  10)  de¬ 
scribes  completely  the  input  for  each 
calculation. 

ORGANIZATION  OF  THE  CALCULATIONS 

Figure  13  shows  the  organization  of  the 
calculations.  The  values  of  all  four  in¬ 
dependent  variables  (p^,  K,  and  failure 
model)  are  shown  for  each  calculation. 

This  was  accomplished  by  constructing  a 
series  of  small  two-dimensional  grids 
(in  pQ  and  S^)  and  then  overlaying  a  large 
two-dimensional  grid  in  K  and  failure  model. 
The  small  grids  within  the  large  grids  are 
positioned  so  as  to  reflect  the  correct 
value  of  K  and  the  correct  failure  model. 

The  locations  of  the  X’s  within  each  small 
grid  give  the  values  of  p^  and  S^.  Thus 
the  calculations  needed  to  investigate  the 
effect  of  a  particular  independent  variable 
can  be  readily  identified.  For  example, 
the  effect  of  density  can  be  determined 
with  the  calculations  performed  for  K  = 

5.0  GPa  and  for  all  failure  models  at  S  = 

w 

0.9.  In  addition,  it  can  be  determined  for 

the  UM-SM  failure  curves  for  S  =0.9, 

w 

0.95,  and  1.0. 

Figure  14  shows  the  calculations  per¬ 
formed  in  (pQ,  S^)  space,  overlaid  on 
Fig.  3.  This  shows  that  the  limits  of  the 


0,7  0.8  0.9  1.0 


Water  saturation 


Fig.  14.  Plot  of  calculations  performed 
in  pQ-S^  space. 

study  are  covered  sufficiently  in  that 
space,  for  all  failure  curves  typical  of 
the  13  Fort  Polk  units.  It  can  also  be 
seen  that  most  of  the  calculations  per¬ 
formed  were  at  S  -  0.9.  Thus  the  study 
w 

was  biased  toward  the  region  where  most 
of  the  Fort  Polk  data  is  concentrated. 


Calculational  Results 


PHENOMENOLOGY  OF  THE  NUCLEAR  DETONATION  to  expand  rapidly,  compressing  the  surround¬ 

ing  soil.  For  a  short  time  after  energy 

The  total  yield,  84  GJ  (20  tons),  re-  release  (^"19  ysec)  the  shock  wave  generated 

sides  initially  in  the  sphere  of  iron  gas  in  the  soil  by  the  iron-gas  expansion  is 

representing  the  nuclear  device.  The  high  strong  enough  to  vaporize  the  soil  (pres- 

initial  pressure  in  the  iron  gas  causes  it  sures  exceeding  110  GPa) .  The  shock  wave 


-17- 


attenuates  very  rapidly,  and  by  the  time 
it  has  reached  a  range  of  approximately 
0.55  m,  it  is  unable  to  cause  further  soil 
vaporization.  Since  this  radius  is 
roughly  twice  the  initial  radius  of  the 
iron-gas  region  (0.27  m) ,  the  total  mass 
of  vaporized  soil  is  considerably  greater 
(roughly  by  a  factor  of  10)  than  the  mass 
of  the  nuclear  source.  The  combined  mass 


of  iron  gas  and  vaporized  soil,  then, 
comprises  the  nuclear  cavity,  which  con¬ 
tinues  to  expand  into  the  surrounding  soil. 

Spherical  divergence  and  inelastic 
effects  continue  to  attenuate  the  peak 
pressures  as  the  shock  wave  travels  be¬ 
yond  the  soil-vaporization  radius.  In 
the  pressure  regions  between  10  and  100 
GPa,  the  shock  wave  can  vaporize  the  free 


Table  8.  Figures  giving  plots  of  kinetic  energy  vs  time  and  peak 
stress  and  particle  velocity  vs  range  for  the  25  SOC74 
calculations. 


Calculation 

Kinetic  energy 

Peak 

stress 

Peak  particle 
velocity 

UM-l.7/0.8 

55 

56 

57 

UM-1.7/0.9 

26. 

45, 

55 

27, 

46, 

56 

28, 

47,  57 

UM-1.7/0.95 

48, 

55 

49, 

56 

50,  57 

SM-1. 7/1.0 

20, 

51, 

55 

21, 

52, 

56 

22, 

53,  57 

UM-1. 9/0.9 

45. 

58 

56, 

59 

47,  60 

A-  UM-1. 9/0. 95 

48, 

58 

49, 

59 

50,  60 

SM-1. 9/1.0 

51, 

58 

52, 

59 

53,  60 

UM-2. 1/0.8 

61 

62 

63 

UM-2.1/0.9 

29, 

45, 

61 

30, 

46, 

62 

31, 

47,  63 

UM-2. 1/0. 95 

48, 

61 

49, 

62 

50,  63 

SM-2. 1/1.0 

23. 

51, 

61 

24, 

52, 

62 

25, 

53,  63 

UX-1.7/0.9 

26 

27 

28 

UX-2.1/0.9 

29 

30 

31 

B- 

SX-1. 7/1.0 

20 

21 

22 

SX-2. 1/1.0 

23 

24 

25 

CM-1.7/0.9 

16, 

32 

17, 

33 

18,  34 

CX-1.7/0.9 

32 

33 

34 

C- 

CM-2.1/0.9 

35 

36 

37 

CX-2.1/0.9 

35 

36 

37 

RM-1.7/0.9 

38 

39 

40 

RX-1.7/0.9 

38 

39 

40 

D- 

RM-2.1/0.9 

41 

42 

43 

RX-2.1/0.9 

41 

42 

43 

CX-1.7/0.9(2.5) 

16 

17 

18 

CX-1.7/0.9(7.5) 

16 

17 

18 
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water  contained  in  the  soil  (but  not  the 
solid  material) .  As  cavity  expansion 
continues,  the  gas  pressure  drops  and  radi¬ 
al  expansion  decelerates.  A  typical  plot 
of  cavity  pressure  versus  time  is  shown  in 
Fig.  15.  By  the  time  the  pressure  at  the 
front  drops  below  10  GPa,  it  has  separated 
from  the  cavity  and  has  started  to  propa¬ 
gate  in  a  manner  similar  to  a  shock  wave 
driven  by  high  explosive  in  the  same  soil. 

EFFECT  OF  INDEPENDENT  VARIABLES  ON  GROUND 
MOTION  AND  ENERGY  COUPLING 

This  section  reviews  the  direct  results 
obtained  from  25  SOC74  calculations.  Many 
quantities  can  be  obtained  from  SOC74,  as  a 
function  either  of  range  or  of  time.  The 
results  presented  in  this  section  are  limit¬ 
ed  mainly  to  kinetic  energy  coupled  to  the 
soil  and  to  peak  stress  and  particle  veloc¬ 
ity  versus  range  in  the  soil.  The  calcu- 
lational  results  are  compared  in  order  to 
show  the  effect  of  each  of  the  four  inde¬ 
pendent  variables:  K,  failure  envelope, 

Pq,  and  S^. 

Location  of  Plots  for  Individual  Calcula¬ 
tions 

In  order  to  allow  for  comparisons  other 
than  those  presented  in  this  section, 

Table  8  lists,  for  each  of  the  25  SOC74 
calculations,  the  figures  that  can  be 
consulted  to  obtain  the  kinetic  energy 
versus  time  and  peak  stress  and  particle 
velocity  versus  range.  They  are  grouped 
according  to  the  letter  designations  ”A” 
through  "E",  as  shown  in  Fig.  13. 

Effect  of  Initial  Bulk  Modulus 

As  already  mentioned,  the  effect  of  the 
initial  bulk  modulus  K  on  ground  motion  was 


Fig.  15.  Typical  curve  of  cavity  pressure 
vs  time  for  an  84-GJ  (20-ton) 
nuclear  detonation,  Z  =  20%. 

felt  to  be  minimal.  Therefore  only  two 
calculations  were  performed  with  values 
that  differed  from  5.0  GPa:  CX-1.7/0.9 
(2.5)  and  CX-1.7/0.9  (7.5),  where  K  was 
equal  to  2,5  and  7,5  GPa,  respectively. 
Unsaturated,  cemented  hypothetical  soils 
were  chosen  for  the  variation  of  K  in 
order  to  maximize  the  effect  for  soils 
typical  of  Fort  Polk.  Figure  16  shows  the 
kinetic  energy  coupled  to  the  soil  for  the 
two  calculations  referred  to  above  and 
CX-1.7/0.9  (this  calculation  according  to 
the  convention  used  in  this  report,  has  K  = 
5.0  GPa).  The  three  curves  are  nearly  the 
same,  indicating  that  varying  K  does 
not  affect  the  coupling  of  kinetic  energy. 
Cratering  is  directly  related  to  the  cou¬ 
pled  kinetic  energy,  and  hence  K  also  does 
not  affect  cratering.  Figures  17  and  18 
show,  for  the  same  three  calculations,  the 
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Fig.  19.  Effect  of  K  on  first  arrival  and  peak  arrival  times. 


peak  stress  versus  range  and  peak  particle 
velocity  versus  range,  respectively. 

Again,  varying  K  has  no  effect  on  the  cal¬ 
culated  results.  Figure  19  shows  that 
varying  K  also  does  not  affect  the  peak 
arrival  time.  First  arrival,  however,  is 
affected  because  K  determines  the  velocity 
of  the  acoustic  wave.  Very  little  kinetic 
energy  is  associated  with  this  "precursor,” 
however,  and  consequently  it  has  very 
little  effect  on  cratering  or  close-in 
ground  motion. 

To  summarize,  for  soils  of  the  Fort 
Polk  type,  K  has  only  a  minor  effect  on 
the  close-in  ground  motion  and  dynamic 
crater  formation  resulting  from  a  nuclear 
detonation. 


Effect  of  Failure 

In  this  section  we  analyze  in  two  ways 
the  effect  of  failure  on  energy  coupling 
and  ground  motion.  First,  we  compare  the 
relative  effect  of  the  maximum  and  mini¬ 
mum  failure  curves  (Figs.  7  through  10)  for 
each  category  (e.g.,  uncemented,  saturated 
soils).  Second,  we  attempt  to  give  the 
total,  qualitative  effect  of  failure  by 
summarizing  results  from  all  of  the  cal¬ 
culations. 

Figures  20,  21,  and  22  show  the  effect 
of  the  maximum  and  minimum  failure  curves 
on  coupled  kinetic  energy  versus  time,  and 
peak  stress  and  peak  particle  velocity 
versus  range,  respectively,  for  the  un¬ 
cemented,  saturated  soils  (SM  and  SX)  with 
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Fig.  20.  Effect  of  failure  on  coupled  kinetic  energy  for  uncemented,  saturated  soils 
(Pq  =  1.7  Mg/m3). 
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Fig.  21.  Effect  of  failure  on  peak  stress 
vs  range  for  uncemented,  satu^ 
rated  soils  (p^  =  1.7  Mg/m^  )  . 


Range  (m) 

Fig.  22.  Effect  of  failure  on  peak  parti*- 
cle  velocity  vs  range  for  unce¬ 
mented,  saturated  soils 
(Pq  =  1. 7  Mg/m3) . 
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Fig.  23.  Effect  of  failure  on  coupled  kinetic  energy  for  uncemented,  saturated  soils 
(Pq  =  2.1  Mg/m^). 


3 

an  initial  bulk  density  =  1.7  Mg/m  , 
Figures  23,  24,  and  25  give  the  same  com- 

3 

parison  for  soils  with  =  2.1  Mg/m  .  In 
all  cases  there  are  only  very  minimal  dif¬ 
ferences  between  the  curves,  indicating 
that,  for  the  saturated  soils,  either 
failure  curve  (or  any  drawn  between  the 
SM  and  SX  curves  in  Fig.  7)  will  adequately 
represent  the  entire  category. 

The  same  comparison  is  shown  in  Figs. 

3 

26,  27,  and  28  for  p^  =  1.7  Mg/m  and  in 
Figs.  29,  30,  and  31  for  p^  =  2.1  Mg/m^ 
for  the  uncemented,  unsaturated  (UM  and 
UX)  category.  These  calculational  re¬ 
sults  show  that,  although  the  differences 

Fig.  24.  Effect  of  failure  on  peak  stress 
vs  range  for  uncemented,  satu¬ 
rated  soils  (p^  =  2.1  Mg/m^) . 
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between  the  peak-stress  and  peak-particle- 
velocity  curves  are  very  small  at  both 
densities,  the  UX  failure  curve  gives  a 
lower  coupled  kinetic  energy  than  does 
the  UM  curve  at  late  times  (>10  msec) . 

The  above  is  seen  to  be  true  also  for 
the  cemented,  unsaturated  soils,  as  shown 

in  Figs.  32,  33,  and  34  for  =  1.7  Mg/ 

3  ^ 

m  and  in  Figs.  35,  36,  and  37  for 

3 

Pq  =  2.1  Mg/m  .  Again,  for  both  densi¬ 
ties,  the  calculated  peak  stresses  and 
particle  velocities  are  not  significantly 
affected  by  the  difference  in  shear 
strength  between  the  CM  and  CX  curves 
(Fig.  9). 

For  the  weak  rocks  there  is  a  large 
difference  between  the  maximum  shear 
strengths  of  the  RM  and  RX  curves  (Fig. 
10).  This  difference,  about  5  MPa,  leads 
to  a  more  substantial  effect  on  kinetic 


Fig.  33.  Effect  of  failure  on  peak  stress 
vs  range  for  unsaturated, 
cemented  soils  (S^  =  0.9, 

Pq  =  1.7  Mg/m^). 


Range  (m) 


Fig.  34.  Effect  of  failure  on  peak  par¬ 
ticle  velocity  vs  range  for 
unsaturated,  cemented  soils 
(S^  =  0.9,  Pq  =  1.7  Mg/m^). 
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Fig,  35.  Effect  of  failure  on  kinetic  energy  coupling  for  unsaturated, 
cemented  soils  (S  =  0.9,  =  2,1  Mg/m^) , 


Fig.  36.  Effect  of  failure  on  peak  stress 
vs  range  for  unsaturated, 
uncemented  soils  (S^  =  0.9, 

=  2.1  Mg/m3). 


Fig.  37.  Effect  of  failure  on  peak  par 
tide  velocity  vs  range  for 
unsaturated,  cemented  soils 
(S„  =  0.9,  p„  =  2.1  Mg/m3). 


30 


Fig.  38.  Effect  of  failure  on  kinetic  energy  coupling  for  weak  rocks 
(S^  =  0.9,  pQ  =  1.7  Mg/m^). 


energy  coupling  and  to  a  slight,  but  no¬ 
ticeable,  effect  on  the  peak  particle  ve¬ 
locity  versus  range.  These  effects  are 

shown  in  Figs.  38,  39,  and  40  for  p  = 

3  ^ 

1.7  Mg/m  and  in  Figs.  41,  42,  and  43  for 

Pq  =  2.1  Mg/m^. 

Thus  increasing  the  strength  within 
any  one  category  lowers  the  kinetic  energy 
coupled  to  the  material  at  late  times 
(>10  msec)  and  also  lowers  the  peak  parti¬ 
cle  velocities  in  the  intermediate  range 
(^10  m) .  The  peak  stress,  however,  is  not 
affected. 

It  is  more  interesting  to  look  at  the 
total  effect  of  failure,  using  all  of  the 
calculations.  In  order  to  do  this,  the 


Fig.  39.  Effect  of  failure  on  peak  stress 
vs  range  for  weak  rocks  ' 

(s„  =  0.9,  Pn  =  1.7  Mg/m3). 
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Fig.  40.  Effect  of  failure  on  peak  par¬ 
ticle  velocity  vs  range  for  weak 
rocks  (S^  =  0.9,  P-  =  1.7  Mg/m^) . 
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Fig.  41.  Effect  of  failure  on  kinetic  energy  coupling  for  weak  rocks  (S  =  0 
p.  =  2.1  Mg/m3) .  ” 


Fig.  42.  Effect  of  failure  on  peak  stress 
vs  range  for  weak  rocks  (S^  =  0.9, 

p-  =  2.1  Mg/m3). 


Range  (m) 

kinetic  energies  were  plotted  at  a  constant 
time  (15  msec)  and  the  peak  particle  veloc¬ 
ities  and  stresses  at  a  constant  range 
(12  m) .  This  time  and  range  were  chosen 
because  they  are  very  important  to  low- 
yield  nuclear  cratering  at  the  optimum 
depth  of  burial  for  soils.  Figure  44a 
plots  the  coupled  kinetic  energy  at  15 
msec  versus  T  .  For  both  sets  of  calcu- 

X  3 

lations  (p^  =1.7  and  =  2.1  Mg/m  )  the 
coupled  kinetic  energy  decreases  with  in¬ 
creasing  T^.  Also,  the  difference  between 
the  two  curves  is  roughly  constant.  Fig¬ 
ures  44b  and  c  plot  the  peak  stress  and 
particle  velocity  at  the  12-m  range  ver¬ 
sus  Although  the  peak  particle  ve¬ 

locity  decreases  slowly  with  increasing 
T  ,  no  clear  correlation  exists  for  the 

X 

two  different  densities,  as  was  seen  for 
the  kinetic  energies  in  Fig.  44a.  The 
peak  stress  seems  to  increase  with  T  , 
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Fig.  43.  Effect  of  failure  on  peak  particle 
velocity  vs  range  for  weak  rocks 
(S  =  0.9,  pf,  =  2.1  Mg/in3)  . 
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Fig.  44,  Total  effect  of  Tx  on  kinetic 

energy  coupling  and  ground  motion. 


but  this  may  not  be  true  in  light  of  the 
artificial  viscosity  Q  used  in  the  SOC74 
code.  The  effect  of  this  term  in  the 
finite-difference  equations  is  to  spread 
the  shock  front  out  in  both  time  and 
space.  In  order  for  the  total  impulse 
to  be  correct,  then,  the  peak  value  must 
be  reduced  somewhat  from  its  "true”  val¬ 
ue.  As  Q  is  increased,  the  effect  in¬ 
creases.  Since  higher  values  of  Q  are 
required  for  the  lower  strength  soils 
than  for  the  weak  rocks  (t^  >  5.0  MPa), 
the  calculated  peak-stress  values  at  the 
lower  values  of  on  both  the  =  1.7 
and  Pq  =  2.1  Mg/m^  curves  are  actually 
slightly  higher.  This  indicates,  then, 
that  there  is  essentially  no  variation  of 
peak  stress  with  T  . 

Effect  of  Density 

As  can  be  seen  from  Fig.  13,  calcula¬ 
tions  were  performed  for  hypothetical 
soils  having  initial  bulk  densities  of 
1.7,  1.9,  and  2.1  Mg/m^.  In  the  pre¬ 
ceding  section  it  was  found  that  both 
peak  stress  and  peak  particle  velocity 
were  not  affected  strongly  by  the  maxi¬ 
mum  shear  strength.  The  kinetic  energy 
was  affected,  but  the  effect  was  similar 

for  both  of  the  densities  (p  =1.7  and 
3  ^ 

pQ  =  2.1  Mg/m  )  used.  In  order  to  study 
the  effect  of  density,  then,  the  failure 
curve,  as  well  as  K,  could  be  held  con¬ 
stant.  For  the  unsaturated  soils,  the 
UM  curve  (Fig,  8)  was  used;  for  the  sat¬ 
urated  soils,  the  SM  curve  (Fig.  7)  was 
used.  By  using  the  results  given  in  the 
preceding  section,  the  effect  of  density 
can  be  found  for  all  failure  models  within 
the  scope  of  the  study. 

Figure  45  shows  the  effect  of  p^  on 
kinetic  energy  coupling  versus  time  for 
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Fig,  45.  Effect  of  density  on  kinetic  energy  coupling  (uncemented  soils,  S  =09) 

w  ^  ^  • 


”  0.9.  Generally,  increasing  in¬ 
creases  the  kinetic  energy  coupled  to  the 
soil.  Figures  46  and  47  show  the  effect 
of  Pq  on  peak  stress  and  peak  particle 
velocity  versus  range,  respectively  (also 
=  0,9).  When  p^  is  increased,  the 
peak  stress  also  increases,  but  the  peak 
particle  velocity  decreases  slightly. 

Figures  48,  49,  and  50  show  the  same 
comparison  for  =  0.95  and  Figs.  51,  52, 
and  53  for  =  1.0,  The  general  effect 
is  the  same  at  all  three  saturations: 
increasing  p^  increases  the  coupled  kinetic 
energy  and  peak  stress,  and  decreases 
slightly  the  peak  particle  velocity. 

In  order  to  examine  the  total  effect  of 
density.  Fig.  54a  shows  the  coupled  kinetic 

Fig.  46.  Effect  of  density  on  peak  stress 

vs  range  (uncemented  soils, 

S  =  0.9). 
w 
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Fig.  52.  Effect  of  density  on  peak  stress  Fig.  53.  Effect  of  density  on  peak  particle 
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Fig.  54.  Total  effect  of  on  kinetic  energy  coupling  and  ground  motion. 


energy  at  15  msec  versus  p^,  and  Figs.  54b 
and  c  show  the  peak  stress  and  peak  parti¬ 
cle  velocity,  respectively,  at  a  range  of 
12  m.  At  all  three  values  of  S^,  the 
kinetic  energy  increases  linearly  with  p^. 
The  coupled  kinetic  energy  at  any  time  is 
the  integral  over  all  space  of  the  density 
p  times  the  square  of  the  particle  veloc¬ 
ity.  By  increasing  p^,  the  quantity 
directly  affected  is  p;  hence  one  might 
expect  a  linear  increase  in  the  coupled 
kinetic  energy  with  increasing  p^.  Figure 
54b  shows  that  the  peak  stress  is  also 
almost  linearly  proportional  to  p^,  espe¬ 
cially  for  =  1.0.  The  peak-particle- 
velocity  plot  (Fig.  54c)  shows  that  the 
effect  of  Pq  is  strong  only  for  ^  1.0. 

In  conclusion,  knowledge  of  the  initial 
bulk  density  p^  is  shown  calculationally 
to  be  important;  it  is  a  convenient  vari¬ 
able  because  all  three  calculational  quan¬ 


tities  examined  are  roughly  linearly 
proportional  to  it. 

Effect  of  Saturation 

The  effects  of  three  of  the  four  inde¬ 
pendent  variables  have  been  discussed  so 
far  in  this  section.  They  have  been  dis¬ 
cussed  in  their  order  of  increasing  impor¬ 
tance  to  ground  motion  and  cratering:  K, 
failure,  and  density.  This  section  re¬ 
views  the  effect  of  saturation  at  three 

different  initial  bulk  densities:  p^  = 

3  ^ 

1.7,  1.9,  and  2.1  Mg/m  .  As  in  the  pre¬ 
ceding  section,  the  comparison  is  made 
for  the  UM  and  SM  failure  models  and  for 
K  =  5.0  GPa. 

Figure  55  shows  the  effect  of  varying 

S  from  0.8  to  1.0  on  kinetic  energy  cou- 
w  3 

pling  for  Pq  =  1.7  Mg/m  .  The  overall 

effect  is  substantial,  especially  at  late 

time  (15  to  20  msec) ,  with  the  kinetic 


Time  (msec) 


Fig,  55.  Effect  of  saturation  on  kinetic  energy  coupling  (uncemented  soils, 
Pq  =  1.7  Mg/m^), 
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Fig.  59.  Effect  of  saturation  on  peak 
stress  vs  range  (uncemented 
soils,  pQ  =  1.9  Mg/m^) , 

energy  coupled  increasing  almost  twofold 

as  S  goes  from  0.8  to  1.0.  Furthermore, 
w 

the  effect  of  S  increases  as  S  ap- 
w  w  ^ 

proaches  1.0:  the  effect  as  goes  from 
0.95  to  1.0  is  roughly  the  same  as  the  ef¬ 
fect  as  S  goes  from  0.8  to  0.95.  Thus 
w 

increasing  increases  the  coupled  kinet¬ 
ic  energy  substantially,  and  it  becomes 
particularly  important  to  know  the  value 

of  S  if  it  is  in  the  range  from  0.95  to 
w 

1.0. 

Figures  56  and  57  show  the  effect  of 

S  on  peak  stress  versus  range  and  peak 
w 

particle  velocity  versus  range,  respec- 

3 

tively,  also  for  =  1.7  Mg/m  .  The 

peak  stress,  like  the  kinetic  energy,  is 

affected  very  strongly  by  S^,  especially 

as  S  goes  from  0.95  to  1.0.  For  example, 
w 


at  a  range  of  10  m,  the  peak  stress  goes 

from  0.012  GPa  at  S  =  0.8  to  0.05  GPa  at 
w 

S  =  0.95,  a  fourfold  increase.  At  S  = 
w  w 

1.0,  the  peak  stress  is  0.17  GPa,  an  in¬ 
crease  of  more  than  threefold  over  its 

value  at  S  =0.95.  Beyond  the  12-  to  15- 
w 

m  range  for  the  unsaturated  soils,  the 
peak  stress  attenuates  very  rapidly. 

This  is  because  the  Fort  Polk  soils  are 
so  compressible  at  low  stresses.  The 
peak-particle-velocity  plot.  Fig.  57,  does 
not  show  as  great  a  variation  as  is  seen 
for  the  peak  stress  out  to  a  range  of 
about  11  m.  Beyond  that  range,  the  un¬ 
saturated  soils  (0.8  <  S  <  0.95)  atten- 
—  w  — 

uate  the  peak  particle  velocity  very 
rapidly. 

Figures  58,  59,  and  60  show  the  kinetic 
energy  coupling  versus  time  and  peak  stress 


Fig.  60.  Effect  of  saturation  on  peak 
particle  velocity  vs  range 
(uncemented  soils, 

Pq  =  1.9  Mg/m3). 
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Fig.  64.  Total  effect  of  on  kinetic  energy  coupling  and  ground  motion. 


and  particle  velocity  versus  range,  respec¬ 
tively,  for  pQ  =  1.9  Mg/m^;  likewise 
Figs.  61,  62,  and  63  show  the  same  calcu- 

3 

lational  results  for  =  2.1  Mg/m  .  The 
general  trends  found  for  the  low-density 
soils  are  followed  also  for  those  at 
higher  densities. 

Figure  64  reviews  the  total  effect  of 

on  kinetic  energy  coupling  and  ground 

motion  for  all  three  densities  considered. 

Thus  the  relative  effect  of  both  p^  and 

S  can  be  seen  in  Fig.  64a  for  coupled 
w 

kinetic  energy.  As  approaches  1.0, 

Pq  becomes  relatively  unimportant.  Over 

:he  relevant  Fort  Polk  range  (0.9  ^ 

S  <  1.0),  S  seems  to  be  the  most  criti- 
w  —  w 

cal  factor  in  the  determination  of  the 


coupled  kinetic  energy.  For  <0.9, 

Pq  and  are  equally  important.  Figure 
64b  shows  the  peak  stress  at  the  12-m  range 
versus  S^.  Again,  in  the  relevant  Fort 
Polk  range  for  S^,  the  quantity  is 
clearly  the  critical  factor  in  determin¬ 
ing  the  peak  stress.  On  the  scale  of  the 
effect  of  on  peak  particle  velocity 
(Fig.  64c)  the  determining  factor  is  only 

S  ,  and  it  is  critical  to  know  its  value 
w 

if  it  lies  between  0.95  and  1.0. 

To  summarize,  then,  S  has  been  found 

to  be  the  most  important  of  the  independent 

bulk  parameters  considered  in  the  study. 

Increasing  S  increases  both  the  strength 
w 

of  the  ground  motion  and  the  kinetic  en¬ 
ergy  coupling,  especially  as  approaches 

1.0. 
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Summary  and  Conclusions 


A  calculational  parameter  study  was 
performed  to  investigate  the  sensitivity 
of  ground  motion  and  energy  coupling  in¬ 
duced  by  a  buried  nuclear  explosive  to 
changes  in  the  bulk  properties  of  the  em¬ 
placement  medium.  Of  specific  interest  in 
the  study  were  soils  typical  of  the  sites 
of  the  ESSEX  I  high-explosive-cratering 
experiments  located  at  Fort  Polk,  Louisi¬ 
ana,  The  calculations  were  one-dimensional 
and  utilized  a  hypothetical  84-GJ  (20-ton) 
nuclear  source.  Constitutive  relations 
based  on  the  bulk  properties  that  were 
varied  in  the  study  were  generated  using  a 
constitutive  model  developed  during  the 
course  of  the  study.  Output  from  this  mod¬ 
el  is  used  directly  as  input  to  SOC74,  the 
computer  code  used  in  the  study. 

The  bulk  properties  varied  included  the 
initial  bulk  density  p^,  the  initial  bulk 
modulus  K,  and  the  initial  water  satura¬ 
tion  S  .  An  additional  property,  failure, 
w 

took  into  account  the  shear  strength  of 
the  soil.  These  properties,  along  with 
the  grain  density  and  Poisson’s  ratio 
(held  constant  throughout  the  study),  made 
up  an  independent  set  from  which  all  other 
bulk  properties  could  be  calculated. 

Results  from  25  SOC74  calculations  were 
plotted,  and  the  relative  effect  of  the 
four  independent  variables  was  examined. 


It  was  found  that,  in  the  order  of  in¬ 
creasing  importance,  K,  failure,  p^,  and 

S  affected  both  the  ground  motion  and  the 
w 

coupling  of  kinetic  energy  to  the  soil. 

Furthermore,  it  was  found  that,  for  values 

of  S  between  0.95  and  1,0,  the  calcula- 
w 

tional  results  varied  considerably,  with 
a  slight  increase  in  leading  to  much 
stronger  ground  motions  and  much  more  ef¬ 
ficient  energy  coupling. 

In  a  study  of  the  cratering  efficiency 
of  explosions  in  various  rock  types  by 
Terhune  et  al.^^  it  was  concluded  that  the 
material  properties  in  order  of  decreasing 
importance  were  water  content,  shear 
strength,  porosity,  and  compressibility. 
Although  the  selection  of  independent  pa¬ 
rameters  in  our  soil  study  varied  somewhat 
from  those  chosen  by  Terhune  et  al.,  the 
results  from  our  study  provide  similar 
conclusions  for  the  prescribed  range  of 
Fort  Polk  soil  parameters.  Water  satura¬ 
tion  is  the  most  sensitive  parameter.  Max¬ 
imum  shear  strength,  which  ranged  from 
about  0.2  to  10  MPa  for  the  soils,  is  of 
less  significance  for  soils  than  it  is 
for  rocks.  For  soils,  bulk  density  (poros¬ 
ity  for  a  given  grain  density)  is  a  more 
important  parameter  to  cratering  than  is 
shear  strength.  Finally,  the  initial  or 
low-pressure  bulk  modulus  is  the  least 
sensitive  parameter  studied. 
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Appendix 


Compressibility  of  Hypothetical  Materials 


Compressibility  was  calculated  for 
each  hypothetical  material  using  the  model 
described  under  "Constitutive  Model  for 
Soils."  The  inputs  for  this  model  were 
the  independent  bulk  properties  plus  the 
merge  and  transition  pressures  associated 
with  the  particular  failure  model  used. 
Values  for  the  independent  bulk  parameters 
and  the  failure  model  used  can  be  obtained 
directly  from  the  name  listed  at  the  top 
of  each  table,  as  discussed  under  "Bulk 
Property  and  Constitutive-Relation  Input 
for  Individual  Calculations,"  while  values 
of  and  P^  are  taken  from  Table  7.  All 


compressibilities  for  hypothetical  ma¬ 
terials  were  derived  assuming  a  mixture 
of  25  wt%  clay  and  75  wt%  sand.  The 
following  tables  list  the  P-y  data  ex¬ 
actly  as  input  to  the  SOC74  code,  for 
all  25  hypothetical  soils  or  weak  rocks 
used  in  the  study.  Since  the  SOC74  pro¬ 
gram  uses  pressures  (P)  in  megabar  units 
and. volumes  (V)  in  cubic  centimeter  units, 
the  actual  tables  are  in  those  units.  The 
quantities  "EF"  and  "EV"  are  the  calcu¬ 
lated  specific  energies  required  to  melt 
and  vaporize  the  hypothetical  soil,  re¬ 
spectively. 
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1.7/0.80  (Fig.  8). 
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Table  A2.  Coordinated  failure  model,  UM  -  1.7/0.90  (Fig.  8) 
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Table  A3.  Coordinated  failure  model,  UM  -  1.7/0.95  (Fig.  8). 
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0.547309 

0 . 480000 

0t .  3  01 2510 

-4.917424 

0.944512 

0.816437 

0.500000 

0.287875 

-6.832386 

1.043368 

1.011575 

0.600003 

0.275745 

-8.243643 

1.133260 

1.112448 

0.300000 

0.256266 

-10.267588 

1.295409 

1.233433 

1 .000000 

0 . 240S43 

-12.967533 

1.442403 

1.360602 

1.500000 

0.211998 

-17.333922 

1.774725 

1.504564 

2.000000 

0. 190775 

-23.559914 

2.083395 

1.619853 

0. 

U .  00000 1 
0 .  ririni'j  ]  n 
0  .  n0OFL.^G 
0 .  i''0n!;,i.::ij 
0 . 00 0000 
U. 000070 
0 ,0001 00 
0.000104 


IJULCiftDING  PATH 

0.071048 
U. 071833 
0.071496 
n .071 344 
0.371 192 
0.370891 
0.970391 
0.970146 
0 . 5 1  ■  0037 


0 


0.029016 
0.029043 
29?90 
0.029564 
0.029337 
0.030381 
0.030922 
0.031728 
0.831835 


0. 

0.036365 

0.036416 

0.036514 

0.036618 

0.036773 

0.036979 

0.037237 

0.037413 
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Table  A4.  Coordinated  failure  model,  SM  -  1. 7/1.0  (Fig.  7) 


Input  variables 


1.7  Mg/m^ 

S  =  1.0 
w 

K  =  5.0  GPa 

10  kPa 

-  311  kPa 

HI 

Output 

0.3417 

(|)q  =  0.5808 

(()  =  0. 
a 

1.3857 

EF  =  6.612  kJ 

EV  =  33.304 

LOADING  PATH 

P 

V 

DP/DV 

MU 

DP/DMU 

0. 

0 .  b  o  o  2  3  5 

0. 

0. 

0. 

0.0OCinc;0 

0. 588234 

-0 . 085000 

0.000002 

0 . 050000 

0.000003 

0.588184 

“0.060052 

0.000087 

0.035322 

0.000010 

0 . 58807  1 

“0.060727 

0.000280 

0.035709 

0. 0000:^-0 

0.587306 

“0.060300 

0 . 000560 

0.035794 

0 . 0L1UU3U 

0.587743 

“0.061106 

0.000838 

0.035894 

0 . 0000-, 0 

0.587417 

-0.061414 

0.001393 

0.036046 

0 . 000070 

0,587094 

“0.061826 

0.001945 

0.036247 

Li.  000  100 

0.586612 

“0.062340 

6.002767 

0.036498 

0 . 000200 

0 . 585M42 

“0.063675 

0.005459 

0.037150 

0 . 00 0400 

0.58204b 

-0.066753 

0.010634 

0.038646 

0.000600 

0.573225 

“0.070908 

0.015555 

0.040640 

0.00  1000 

0.574038 

“0.077113 

0.024732 

0.043588 

0.002000 

0.563115 

“0.091550 

0.044609 

0.050309 

0.004000 

0.546851 

“0. 122969 

0.075678 

0.064374 

0.006000 

0.534893 

-“0. 167259 

0.039724 

0.083171 

0.003000 

0.521413 

“0.222552 

0.128155 

0. 105519 

4.00961  ;0 

“n.264r!54 

0.  133078 

0. 121879 

4.0 innrn 

I-.M  ,'744 

“0.276209 

0, 136248 

0.  126199 

o.nu  !>i  n 

IK ' 

“0.3341 13 

0.  170070 

0. 147631 

(4  !  i.'i'i  •!  ii 

4.  141444 

“0.447001 

0.  136696 

0.  187786 

0  •  '  L  M '  M  •  4 

4 . 

“0.546760 

0.214776 

0.221242 

070444 

0.4, '',.471 

“0.652023 

0.238308 

0.254970 

4 . 04I  )44!:1 

0.4t>2713 

“0.811816 

0.271274 

0.303348 

li .  064! 'i  iO 

0.444252 

“1 .083348 

0.324102 

0.378581 

4..  14440!,1 

11.4  0:42  1 

“1.548510 

0 . 405846 

0.439335 

It.  iMtriiio 

iJ.. 4  147  16 

“1.660881 

0.514835 

0.458761 

!  1 .  OltOl  |40 

0.  4/344,' 

-2.221264 

0.603050 

0.536397 

4 .  Olji  i4i:.0 

0.724744 

“2 . 7 463 1 0 

0.735810 

0.562556 

0 . 400000 

0.71,18441 

“4.772514 

0.907125 

0.824237 

0 . 000000 

0.247214 

“6.567338 

1.006164 

1.009706 

0 . 600000 

0.280559 

“7.902333 

1.036651 

1.  105132 

0 . 800008 

0.260176 

“9.811885 

1.260913 

1.217569 

1 .0010000 

0.243382 

“12.350573 

1.410976 

1.332767 

3 .  rjorjnoo 

0.213573 

“16.442434 

1.754258 

1,456531 

2 . 000000 

0.  191106 

“22.254610 

2.078061 

1.544150 

UNLOADING  PATH 

0. 

0.500235 

0.000000 

0. 

0.0O0G00 

0.588234 

0.000003 

0.035643 

0.000003 

G. 588104 

0.000087 

0.035659 
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Table  A5.  Coordinated 

failure  model,  UM  - 

1.9/0.90  (Fig.  8). 

Input  variables 

Pq  =  1.90  Mg/m^ 

S  =  0.90 
w 

K  =  5.0  GPa 

=  100  kPa 

P,,  =  10.4  MPa 

M, 

Output 

Z  =  0.2061 

(i)Q  =  0.4350 

(|)  =  .0435 

a 

e  =  0.6930 
w 

EF  =  7.297  kJ 

EV  =  35.380  kJ 

LOADING  PATH 

F 

V 

DP/DV 

hU 

DP/DMU 

0. 

0.526516 

0. 

0. 

0. 

0.000001 

0.526305 

-0.095002 

0.000020 

0 . 050000 

G. 0000 10 

0.514605 

-0.000769 

0.022756 

0.000396 

0.000020 

0.511108 

-0.002860 

0.029754 

0.001429 

0.000030 

0.509041 

-0.004838 

0.033935 

0.002392 

0.000050 

0.506385 

-0.007530 

0.039358 

0.003688 

0.000070 

0.504582 

-0.011090 

0.043073 

0.005384 

0.000100 

0.502597 

-0.01511 1 

0.047193 

0.007281 

0.000104 

0.502372 

-0.017836 

0.047661 

0.008556 

0.000200 

0.501436 

-0. 102546 

0.049616 

0.049031 

0 . 000400 

0,499572 

-0. 107294 

0.053533 

0.051067 

0 . 000600 

0.497814 

-0. 113750 

0.057254 

0.053749 

0.00 1000 

0.494571 

-0. 123361 

0.064185 

0.057707 

0.002000 

0.487700 

-0. 145520 

0.079130 

0.066689 

0.004000 

0.477326 

-0.192806 

0.102633 

0.035279 

0.006000 

0.469557 

-0.257403 

0.  120878 

0. 109615 

0.009000 

0.460600 

-0.334940 

0. 142675 

0.137636 

0.009600 

0.459069 

-0.392105 

0. 146484 

0.  157528 

0.010000 

0.458087 

-0.407155 

6.143943 

0. 162682 

0.015000 

8.447717 

-0.482135 

0.175556 

0.187877 

0.020000 

0.439656 

-0.620330 

0. 197107 

0.232003 

0.024000 

0.434214 

-0.734986 

0.212112 

0.266593 

0.030000 

0.427152 

-0.849671 

0.232150 

0,299427 

0.040000 

0.417287 

-1.013618 

0.261281 

0.343278 

0.060000 

0.401556 

-1.271413 

0.310690 

0.404782 

0. 100000 

0.377692 

-1.676137 

0.393506 

0.483000 

0,  150000 

0.350947 

-1.869527 

0.439701 

0.470831 

0.200000 

0.330130 

-2.401855 

0.594269 

0.528721 

0.300000 

0.295912 

-2.922465 

0.778621 

0.542439 

0.400000 

0.278116 

-5.619257 

0,892431 

0.878662 

0.500000 

0.265915 

-8.  196100 

0.979261 

1 .  151677 

0,608000 

0.255940 

-10.024623 

1.056404 

1.296292 

0.SO0000 

0.240182 

-12.692395 

1.  191317 

1.482437 

1.000000 

0.227930 

-16.322694 

1.309116 

1.697803 

1.500000 

0.205518 

-22.309494 

1.560929 

1.985605 

2.806000 

0,  189415 

-31.051538 

1.778634 

2.296681 

0. 

UNLOADING  PATH 

0.503420 

0.045481 

0. 

0.000001 

0.503403 

0.045503 

0.047030 

0.000010 

a.  5033 17 

0.045694 

0.047098 

0.006020 

0 . 5032 15 

0.045905 

0,047227 

0.800030 

0,503114 

0.046117 

0.047363 

0 .  eocitefc 

0.502312 

0.046537 

0.047566 

0 . 0060/^3 

0.502711 

0.046955 

0.047838 

0.000100 

0.502412 

0.047578 

0.048177 

0.000104 

0.502372 

0.047661 

0.048407 
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Table  A6.  Coordinated  failure  model,  UM  -  1.9/0.95  (Fig.  8) 


Input  variables 


)q  =  1.9  Mg/m^ 

=  100  kPa 

S  =0.95 
w 

=  10.4  MPa 

M 

Output 

K  = 

=  5.0  GPa 

Z  =  0.2238 

(j)Q  =  0.4477 

=  0.0224 

:  =  0.7700 

w 

EF  =  7.485  kJ 

EV 

=  36.063  kJ 

i 


LOADING  PATH 


P 

n. 

D . Dun no  1 

0.  DO  nil  10 

o. ooDnno 

n.onoono 

D  .  LID  no  no 
0.000070 
0.000100 
0.000104 

0 .  oooznuo 

0 . 000400 
0 . 000600 
0.001000 
0.082000 
0.004000 
0 . 005000 
0 . 009000 
0 . 009600 
0 . 0 1  no  nr 
0 . 0  1  bOOL 
0 . 020000 
0 . 024000 
■0 . 030000 
0 . 040000 
U . 060000 
D.  100000 
0.  130000 
0 . 200008 
G . 30 00 GO 
0 . 400000 
0 . 500000 
0 . 600000 
0 . 800000 
1 . 00 0000 
1.500000 
2.000000 


11.000:10  1 
0 .01'.! ;  1 0 
0 . 00i.:i  ..,'0 
r.i .  ij00u..;o 
U . 000000 
0.000070 
0.000100 
0.ODO1O4 


o.okoiG 

0.526305 
n . 520297 
0.518434 
0.517304 
0.515807 
0.514752 
0.513544 
0.513404 
0 . 512393 
0.510381 
0 . 508483 
0 . 504986 
0.497583 
0.486438 
0.478119 
0 . 468569 
0.466943 
0.465399 
G. 454927 
0.446457 
0.440769 
0.433427 
0.423240 
0.407152 
0.3330.30 
0.355844 
0.334805 
G . 300300 
0 . 282090 
0.263493 
0.253166 
0.242803 
0,230033 
0.206573 
0.  189637 


DP/DV 

0. 

“0.095002 
-0.001493 
-0.005368 
-0.008853 
-0.013353 
-0.013951 
-0.024838 
-0.028593 
-0.094962 
-0.099385 
-0. 105404 
-0. 114371 
“0. 135084 
-0. 179450 
-0.240422 
-0.314149 
-0.368878 
-0.383338 
-0.455679 
-0.590332 
-0,703258 
-0.817207 
-0.981639 
-1.243131 
-1.658226 
“1.839247 
-2.376533 
“2.893103 
-5.491570 
-7 . 933094 
-9.633792 
-12.222165 
-15.662501 
-21.312847 
-29.522196 


MU 

0. 

0.000020 
0.011569 
0.015204 
0.017420 
0.020373 
0-.  022466 
0.024870 
0.025149 
0.027172 
0.031222 
0.035070 
0.042239 
0.057745 
0.031980 
0. 100305 
0. 123240 
0. 127152 
0.  129677 
0.  156924 
0.  178872 
0. 194085 
0.214312 
0,243540 
0.292677 
0.374087 
0.479061 
0.572005 
0.752633 
0.865770 
0.952336 
■1.030303 
1.  167669 
1.287998 
1.547841 
1.775388 


0.022896 
0,0229 18 
0 . 0231  I b 
0.023335 
0.023553 
0 . 023988 
0.024420 
9.025064 
0.025149 


DP/DMU 

0. 

0 . 050000 
0.000779 
0.002751 
0.004511 
0.006772 
0.009560 
0.012475 
0.014326 
0.047464 
0.049382 
0.051974 
0.055799 
0.064491 
0.082526 
0.  106241 
0.  133721 
0. 153347 
0. 158450 
0.  183504 
0.227809 
0.262941 
0.296628 
0.342144 
0.407018 
0.491344 
0.476306 
0.537960 
0 . 553625 
0.883881 
1.146583 
1.285065 
1.461281 
1.662106 
1.924243 
2, 197349 


0. 

0.045484 
0.045549 
0.045674 
0.045805 
0.046002 
0 . 046264 
0.046591 
0.046813 


iiMi  nnniNG  oath 


Mnnnn 

Minnn 

i4.vj  4 
M4n05 
M398G 
5137G9 
5 1 3447 
513404 
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Table  A7.  Coordinated  failure  model,  SM  -  1. 9/1.0  (Fig.  7) 


Input  variables 


pQ  =  1.9  Mg/m^ 
=  10  kPa 

Z  =  0.2427 

e  =  0.8556 
w 


S  «  1.0 
w 

=  311  kPa 

M 

Output 
(|)q  =  0.4611 
EF  =  7.706  kJ 


K  =  5.0  GPa 


(f)  =  0. 
a 

EV  =  36.832  kJ 


LOADING  PATH 

p 

V 

DP/DV 

MU 

DP/DflU 

0. 

0.526316 

0. 

0. 

0. 

0.000000 

0.526315 

-0.095000 

0.080002 

0.050000 

0.080003 

0.526279 

-0.083533 

0.000070 

0.043962 

0.000010 

0 . 526197 

-0.083999 

0.000226 

0.044197 

0.000070 

0.526078 

-0.084234 

0.000452 

0.044304 

0 . 0000:^^0 

0.525960 

-0.034514 

0.600677 

0.044431 

0.000000 

0.525724 

-0.084933 

0.001125 

0.044621 

0.000070 

0 . 525490 

“0.035492 

0.001571 

0.044875 

0 . 000  1. 00 

0.525142 

-0 . 036191 

0.002235 

0.045191 

0.000200. 

0.524006 

-0.083002 

0.004408 

0.046011 

0.000000 

0.521836 

“0.892184 

0.008534 

0.047894 

0.000600 

0.519731 

-0.097799 

0.012552 

0.050402 

0.001000 

0.516024 

-0. 106172 

0.019945 

0.054108 

0.002000 

0.503058 

-0. 125541 

0.035936 

0.062535 

0.004000 

0.496095 

-0. 167176 

0.060918 

0.080858 

0,006000 

0.487194 

-8.224710 

0.080299 

0.  10319.1 

0.0O'=’OO0 

0.477017 

-0.294757 

0. 103349 

0.130153 

0.009600 

0.475288 

-0 . 347087 

0.  107362 

0. 149514 

0.010000 

0.474180 

-0.360959 

0.  109950 

0. 154565 

0.01 5000 

0.462569 

-0.430635 

0.  137811 

0. 179466 

0.020800 

0,453665 

-0.561549 

0. 160142 

0.223899 

0 . 024000 

0.447717 

“0.6724S9 

0.  175555 

0.259524 

0,030000 

0.440078 

-0,785401 

0.  195962 

0.294021 

0.040000 

0.429550 

-0.949877 

0.225273 

0.341 164 

0.060000 

0.413082 

-1.214498 

0.274119 

0.409450 

0.  100000 

0.338687 

-1.639656 

0.354087 

0.500199 

0.  150000 

0,361035 

-1.808206 

0.457797 

0.482115 

0 . 200000 

0 . 33976 1 

“2.350272 

0.549077 

0.547765 

0.300000 

0.304951 

-2,872722 

0.725904 

0.565523 

0.400000 

0.286302 

-5.362423 

0.838321 

0.889548 

0.500000 

0.273285 

-7.681798 

0.925889 

1. 141976 

0 . 600000 

0,262586 

-9.346371 

1.004356 

1.274416 

0.880060 

0.245580 

-11.760374 

1-.  143156 

1.440917 

1.000000 

0,232263 

-15.018677 

1.266034 

1.627638 

1.500000 

0.207692 

“20.349332 

1.534115 

1.865108 

2.000000 

0.  133872 

-28.057565 

1.771956 

2. 102247 

UNLOADING  PATH  . 

0. 

0.526316 

8. 

0. 

0 . 000000 

0.526315 

0.000002 

0.0441 14 

0.000003 

8.526279 

0.000070 

0.044134 
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Table  A8.  Coordinated  failure  model,  UM  -  2.1/0*80  (Fig.  8) 


Pq  =  2.1  Mg/m^ 

=  100  kPa 

Z  =  0,1161 
e  =  0.3508 


Input  variables 

S  =  0.80 
w 

«  10.4  MPa 

M 

Output 
(|)q  =  0.3048 
EF  =  8.226  kJ 


K  =  5.0  GPa 


c()  =  0.0610 

EV  =  37.99  kJ 


< 


LOADING  PATH 


p 

V 

DP/DV 

MU 

DP/DMU 

0. 

0.476190 

0. 

0. 

0. 

0 . 06000 1 

0. 476131 

-0. 105002 

0 , 000020 

0 , 050000 

0.000618 

0.461273 

“0.000604 

0.032326 

0.000279 

0.000020 

0.456927 

“0.002298 

0.042158 

0.001017 

0.000ij30 

0.454394 

“0.003947 

0.047969 

0.001721 

0.000050 

0.451196 

“0 . 006254 

0.055396 

0.002693 

0.000O70 

0.449073 

“0.009420 

0.060386 

0.004008 

0 . 000100 

0.446791 

“0.013150 

0.065800 

0.805541 

0.000104 

0.446533 

“0.015773 

0.066406 

0.006608 

0.000200 

0.4459S0 

“0. 172121 

0.067739 

0.071983 

0.000^100 

0 . 444367 

“0. 17964S 

0.070411 

0.074849 

0 . 1  lOifJbOiO 

0 . 4433 1 3 

“0,  189825 

0.072952 

0.078705 

0 . 00 1000 

0.441861 

“0.204852 

0.077694 

0.084362 

0.002000 

0.437676 

“0.238958 

0.087998 

0.097046 

0 .0040:  . 10 

0,431210 

“0.309336 

0.104311 

0. 122600 

O.O060’:iO 

0.426213 

“0.400621 

0. 117246 

0.  154623 

0.009000 

0.420261 

“0.503644 

0.  133081 

0.  189450 

O.009611U 

0.419219 

“8.575516 

0.  135899 

0.212930 

0.010000 

0.413545 

“0.593924 

0,137727 

0.218843 

0.015000 

0.41 1222 

“0,682765 

0. 157988 

0.246780 

0.020000 

0 . 405235 

“0.835130 

0.  175097 

0.292252 

0. 024000 

0.401036 

“0.952433 

0. 187402 

0.325062 

0.030000 

0.395394 

"1.063508 

0.204344 

0.354138 

0 . 040000 

0.337155 

“1.213749 

0.229974 

0.390178 

0 . 060000 

0.373236 

“1.436364 

0.275844 

0.436017 

0.  100000 

0 . 350676 

“1.773069 

0.357921 

0.487343 

G.  150000 

0.326160 

“2.039471 

0.453991 

0.489862 

0.200000 

0 . 306465 

“2.538763 

0.553816 

0.532909 

0 . 300000 

0.273703 

-3.852330 

8.739805 

0.537665 

0 . 400000 

0.258002 

“6.368760 

0.845687 

0 . 944446 

0,500000 

0.247303 

“9.809391 

0.921611 

1.317108 

0.600000 

0.239610 

“12. 197430 

0.987360 

1.520930 

0 . S0OO00 

0.226921 

“15,761653 

1.098490 

1.799695 

1.006000 

0.217232 

-20.749450 

1.  191581 

2.148438 

1 .500000 

0.280174 

“29.227079 

1.378879 

2.669538 

2.006000 

0.  1S3294 

“42.086405 

1.528974 

3.331239 

UNLOAD FNG  PATH 

0. 

0.447161 

0.064920 

0. 

0.000001 

U. 0000  10 

0.447155 
e. 447 100 

0 . 064335 
0.065065 

8.069013 

0.069111 

0.009020 

LI .  447039 

0.065203 

0,069298 

0.000930 

0.446979 

0 . 065353 

0.069495 

0 . 000050 

0.446853 

0 . 065640 

0.069790 

0.800070 

0.446733 

0.065325 

0.070184 

G . 000100 

0.446561 

0.066343 

0.070674 

0 . 000 104 

0.446538 

0.066406 

0.071008 
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Table  A9.  Coordinated  failure  model,  UM  -  2.1/0.90  (Fig.  8) 


Input  variables 


2.1  Mg/m^ 

S  =  0.90 
w 

K  > 

“  5.0  GPa 

>  - 
T 

100  kPa 

=  10.4  MPa 

M 

Output 

Z  « 

0.1380 

(pQ  =  0.3220 

=  0.0322 

'w 

0.4275 

EF  =  8.551  kJ 

EV 

=  39.038  kJ 

LORDING  FDTH 

P 

DP/DV 

MU 

DP/DMU 

P  . 

o.ir'biso 

0. 

0. 

0. 

u.ODnnoi 

0.476181 

”0. 105002 

0,000020 

0 . 050000 

i.M.u:u:ini0 

U.. 168390 

“0.001155 

0.016653 

0.000541 

i).  nuon.'^o 

U .  •466042 

“0.004253 

0.021777 

0.001952 

0 . 000030 

U . 464643 

“0.007181 

0.024839 

0.003265 

n . nnoDoo 

0.462854 

“0.011 142 

0.028814 

0.005032 

0 .  or;OMr0 

0.461631 

“0.016360 

0.031538 

0.007341 

0.000100 

0 .  ■160282 

“0.022233 

0.034562 

0.009921 

0 .  orio  1 04 

0.460129 

“G. 026201 

0.034905 

0.011653 

[} .  000200 

0 . 453480 

“0. 147723 

0.036369 

0.065586 

0 . 000400 

0.458 134 

“0.154316 

0.039301 

0.068224 

0 . 00 0000 

0.456958 

-0.163243 

0.042087 

0.071777 

0.001000 

0.454692 

“0. 176474 

0.047282 

0.077001 

0.002000 

0.449S53 

“0.206657 

0.058547 

0.088768 

0 . 004000 

0.442436 

“0.269662 

0.076292 

0. 112710 

O.OOGOOO 

0.436768 

“0.352833 

0.090260 

0. 143182 

0.003000 

0.430081 

“0.448637 

0.  107212 

0.  176976 

O.O03GO0 

0.428920 

“0.516684 

0.  1 10209 

0.200157 

0.O10OOO 

0.428171 

“0.534266 

0.  1 12151 

0.206049 

0.015000 

0.420106 

“0.619965 

0.  133501 

0.234187 

0.020000 

0.413614 

“0.770203 

0.  151291 

0.281049 

0.024000 

0.409112 

“0.888482 

0.  163961 

0.315722 

0.030000 

0.403125 

-1.0021 15 

0. 181248 

0.347071 

G . 040000 

0.394490 

“1 .  158088 

0.207105 

0.386755 

0 . 000000 

0.330130 

“1.392745 

0.252705 

0.438590 

0. 100000 

0.357252 

-1.748455 

0.332925 

0.498633 

0.  150000 

0.332194 

“1.995318 

0.433472 

0.497277 

0.200000 

0.312226 

“2.504019 

0.525148 

0.545402 

0 . 300000 

0.279110 

“3.019666 

0.706106 

0.552614 

0.400000 

0.262898 

“6.  168479 

0.811312 

0.950516 

0:500000 

6.252216 

“9.361209 

0.888028 

1.303500 

0.600000 

0.243585 

-11.586166 

0.954927 

1.494794 

0.800000 

0.230149 

-14.885425 

1.069055 

1.752425 

1.000000 

0.219374 

-19.464453 

1.165746 

2.068443 

1.500000 

0.201475 

-27. 175850 

1.363522 

2.528115 

2.000000 

0.  188567 

“38.735521 

1.525313 

3.090405 

UNLOADING  PATH 

0. 

0 . 460R56 

0.033275 

0. 

0.000001 

0.460848 

0.033291 

0.062858 

G. 000010 

0 . 460785 

0.033434 

0.062948 

0.0000^0 

0.460714 

0.033592 

0.063120 

G.0L1GO30 

0 .460644 

0.033750 

0.663301 

0 . 000050 

0 . 460503 

0.034065 

0.063572 

0.000070 

0 . 4b0364 

0.034378 

0.063933 

0 . 000 100 

0.460157 

0.034844 

0.064384 

8.000104 

0.460129 

0 . 034905 

0.064691 
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Table  AlO.  Coordinated  failure  model,  UM  -  2*1/0.95  (Fig.  8) 


Pq  =  2.1  Mg/m^ 

=  100  kPa 

Z  =  0.1499 
e  =  0.4709 


Input  variables 

S  =  0.95 
w 

=  10-4  MPa 

Output 
(f)Q  «  0.3314 
EF  =  8.749  kJ 


K  =  5.0  GPa 


<f)  =  0.0166 

EV  =  39.645  kJ 


i 


LOADING  PATH 


p 

V 

DP/DV 

MU 

DP/DMU 

0. 

0.476190 

0. 

0. 

0. 

0.000001 

0.476181 

-0. 105002 

0.000020 

0.050000 

0.000010 

0.472169 

“0.002243 

0.008518 

0.001059 

0.000020 

0.470318 

“0.007993 

0.011197 

0.003732 

0.000030 

0.470157 

“0.013150 

0.012833 

0.006114 

0.000050 

G. 469 147 

“0.019793 

0.015014 

0.009168 

0 . 000070 

0.468432 

“0.028002 

0.016562 

0.012923 

0.000100 

0.467613 

“0.036608 

0.018343 

0.016839 

0.000104 

0.467518 

“0.042084 

0.018550 

0.619321 

0 . 000200 

0. 4663 18 

“0. 137155 

0.020077 

0.062860 

0 . 000^100 

0,465423 

“0. 143330 

0.023136 

0.065396 

O.OOOiSGO 

0.464104 

“0. 151703 

0.026042 

0.068814 

0.001000 

0.461667 

“0. 1641 15 

0.031459 

0.073844* 

0.002000 

0 . 456472 

“0. 192511 

0.043197 

0.035196 

0 . 004000 

0 . 448533 

-0.252087 

0.061649 

0. 108389 

0 . 0060 00 

0.442503 

“0.331347 

0.076130 

0. 138108 

0.003000 

0.435419 

-0.423493 

0.093638 

0.  171352 

0.003600 

;0. 434193 

“0.489483 

0.096726 

0. 194333 

0.010000 

0 . 433403 

-0.506604 

0.098724 

0.200199 

0.015000 

0.424935 

“0.590442 

0. 120619 

0.228356 

0.020000 

0.418163 

“0.738981 

0.  138751 

0.275757 

0.024000 

0.413503 

“0.857172 

0.  151602 

0.311255 

0.030000 

0.407327 

“0.971625 

0.  169061 

0.343669 

0.040000 

0.398477 

“1.129921 

0.195026 

0.385136 

0.060000 

0.385377 

“1.369880 

0.240476 

0 . 440046 

0.  100000 

0.360827 

“1.735360 

0,319718 

0.504778 

0.  150000 

0.335474 

“1.972110 

0.419457 

0.501312 

0 . 200000 

0.315357 

-2.485527 

0.510003 

0.552204 

0 . 300000 

0.282048 

“3,002201 

0.683329 

0.560771 

.'0.400000 

0.265560 

“6.064802 

0.793157 

0.953942 

■0.500000 

0.254612 

-9.134106 

0.870260 

1.296959 

0.6GOG0O 

0.245746 

“11,273876 

0.937737 

1.482005 

0.800000 

0.231904 

“14.448779 

1.053397 

1.729196 

1.000000 

0.221283 

“18.830524 

1.151956 

2.029256 

1.500000 

0.202132 

“26. 177152 

1.355256 

2.459414 

2.000000 

0.  188715 

“37. 128547 

1.523327 

2.974931 

0. 

UNLOADING  PATH 

0 . 468300 

0.016849 

0. 

0.000001 

0.468292 

0.016866 

0.060238 

0.000010 

0.468224 

0.017815 

0.060325 

0 . 00O0jt^ 

0 . 600OTO 

0.463148 

0.017180 

0.060490 

8.468072 

0.017345 

0.060664 

0 . O0O05G 

0.467921 

0.017673 

0.060924 

0 . 

0,467771 

0.018000 

0.061271 

0.000100 

0.467547 

6.018486 

0.061705 

0.000104 

0.467518 

0.018550 

0.061999 

-56- 


Table  All.  Coordinated  failure  model,  SM  -  2. 1/1.0  (Fig.  7) 


Input  variables 


2.1  Mg/m^ 

10  kPa 

S  «=  1.0 
w 

=  311  kPa 

Output 

K  = 

5.0  GPa 

0.1625 

(j)Q  =  0.3413 

,  = 

0. 

0.5182 

EF  =  8.926  kJ 

EV  = 

40.318 

i  ni  ilUNG  POTH 

!• 

V 

DP/DV 

MU 

DP2DMU 

i1  , 

0.47i;iOf.) 

0. 

0. 

0. 

M-.  1  ii  M'-iuin 

n.47i>i2n 

-0.  105000 

0.000002 

0.050000 

i  1 .  l-!i'iiih  •: 

!).47I>  16!3 

"0.  121532 

0.000854 

0.05?'893 

ij.nrMiUi.t 

11.476100 

-0,  121777 

0.000173 

0.057976 

IJ  .  UL.ll.H 

0.4 76026 

-0.  122109 

0.000345 

0.058117 

I'l .  unou^n 

0 . 475045 

-0. 122503 

0.000516 

0.0532S4 

m . 000050 

0.475732 

-0. 123092 

0.000853 

0.058535 

0.000070 

0.475621 

-0. 123879 

0.001198 

0.058869 

O.onnioo 

0.475330 

-0.  124862 

0.001704 

0.059286 

0 . 000200 

0.474536 

-0. 127403 

0.003361 

0.060365 

0.000400 

0 . 4  7  3  U  3  5 

-0. 133273 

0.006543 

0.062840 

0 . 000600 

0.471673 

-0. 141 124 

0.009567 

0.066132 

0.001000 

0 .  .-IbSOSS 

-0. 152776 

0.015203 

0 . 070982 

0.00.2000 

6.463433 

-0. 179489 

0.027406 

0.081945 

0.004000 

0.455006 

-0.235793 

0,046558 

0. 104423 

0.006000 

0 . 44853 1 

-0.311253 

0.061549 

0. 133413 

0 . 002000 

0.441076 

“0.599743 

0.079611  ' 

0. 166096 

0 . 000600 

0.433732 

“0.463615 

0.082788 

0. 188855 

0.010000 

0.438349 

-0.480250 

0.084843 

0. 194687 

0.015000 

0.430053 

“0.562074 

0. 107282 

0.222818 

0.020000 

0.422396 

-0.708533 

0.  125755 

0.270669 

0.024800 

0.413156 

-0.826309 

0.  138787 

0.306928 

0 . 030000 

0.411 78 1 

-0.941272 

0. 156416 

0.340361 

0 . 040000 

0.402703 

-1  .  101526 

0.  182485 

0.383588 

0 . 060000 

0.337849 

-1.346453 

0.227772 

0.441630 

0 .  100000 

0.364616 

-1.721694 

0.306005 

0.511298 

0 . 150000 

0 . 333950 

-1.948095 

0.404833 

0.505593 

0 . 200000 

0.313676 

“2.466225 

0.494277 

0.559419 

0 . 300000 

0.235163 

-2.383910 

0.669838 

0.569440 

0.400600 

0.263381 

-5.358647 

0.774309 

0.957661 

0.500000 

0.257151 

-3.305133 

0.851791 

1.290625 

0.600000 

0.243036 

-10.970493 

0.319844 

1.469430 

0.800000 

0.233764 

-14.013109 

1.037060 

1.706261 

1.000008 

0.222776 

-18.202203 

1.  137531 

1.990622 

1.500000 

0.202931 

-25. 195779 

1.346559 

2.392019 

2.000000 

0. 138873 

-35.564769 

1.521227 

2.862578 

UHLCADING  PAT 

0. 

0 . 000000 
0.008003 

0.476190 

0 . 476 1 90 

0.476165 

0.000000 

0.000082 

0 . 000054 

0. 

0.057867 

0.057893 
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Table  A12.  Coordinated  failure  model,  UX  -  1.7/0.90  (Fig.  8) 


Input  variables 


‘o  * 

1.7  Mg/m^ 

S  =  0.90 
w 

K  = 

=  5.0  GPa 

T  “ 

100  kPa 

=  14.7  MPa 
.  M 

Output 

Z  = 

0.2901 

(|)q  =*  0.5480 

=  0.0548 

'w 

1.0913 

EF  =  5.938  kJ 

EV 

=  31.773  kJ 

LOADING  PATH 

P 

V 

DP/DV 

MU 

DP/DMU 

0. 

0.588235 

0. 

0. 

0. 

0.000001 

0.538224 

-0.085002 

0.000020 

0 . 050000 

0.000010 

0.572764 

-0.000582 

0.027012 

0.000333 

0.000020 

0.568163 

-0.002173 

0.035329 

0.001202 

0.000030 

0 . 565445 

-0.003680 

0.040304 

0.002010 

0.000050 

0 . 561953 

-0.005727 

0.046769 

0.003094 

0 . 000070 

0.559581 

-0.008430 

0.051207 

0.004506 

0.000100 

0.5569r56 

-0.011473 

0.056142 

0.006079 

0.00014? 

0,553972 

-0.015700 

0.061850 

0.008235 

0.000200 

0.553265 

-0.074946 

0.063207 

0.039050 

0.000400 

0 . 550699 

-0.077952 

0.068160 

0.040376 

0.000600 

0.543233 

-0.032755 

0.072869 

0.042478 

0.001000 

0.543B35 

-0.089927 

0.081644 

0.045584 

0.002000 

0 . 53445 1 

-0.  106573 

0, 100634 

0.052659 

0.004000 

0.520426 

-0. 142600 

0. 130295 

0.067427 

0 . 006000 

0.510060 

-0. 192940 

0.  153266 

0.087067 

0.003000 

0 . 493299 

-0.255634 

0. 180486 

0.  1 10216 

0.009600 

0.496313 

-0.302120 

0.185209 

0.  127021 

0.010000 

0 . 495042 

-0.314674 

0.188253 

0.131434 

0.015000 

0.431824 

-0.378255 

0.220852 

0.  153378 

0.020000 

0.471326 

-0.500113 

0.246721 

0.  193280 

0.024000 

0.465222 

-0.605718 

0.264418 

0.226028 

0.030000 

0.456834 

-0.715261 

0.287636 

0.258424 

0.040000 

0.445448 

-0.878275 

0.320548 

0.303833 

0 . 060000 

0.428024 

-1 . 147882 

0.374303 

0.372059 

0,  100000 

0.402941 

-1.594661 

0 . 459856 

0 . 467549 

0.  150000 

0.374113 

-1.734453 

0.572346 

0 . 444483 

0.200000 

0.352247 

-2.286603 

0.669953 

0.512260 

0 . 

0.:ili;668 

-2.810702 

0.857575 

0.532985 

0 .  i'ij 

n..".H4Ti5 

-5.062453 

0.981156 

0.809183 

0  .!-4 11^0:10 

A '82838 

-7.  103901 

1.079758 

1.014183 

.  hi'i  ■O:  iO 

8.:',M;'il2 

'-8.593845 

1.  168992 

1.  120644 

.0.  :)n; -nri!! 

8 .  '5  :'  u‘r 

-  1  n . 758 1 34 

1.328936 

1.250444 

1  , 

u,.‘ 

“13. 009 1 75 

1.472814 

1.390061 

1 . 

0.2 in:; 1 1 

“18.26C«455 

1.794316 

1.555199 

2 .  ui.ujnoo 

0.  j;jU463 

-24,940076 

2.088445 

1.699939 

UNLOADING  PATH 

0. 

0.1^^^8000 

0 . 05?980 

0. 

0.000001 

0.555085 

0.058007 

0.037216 

0.000010 

0.555858 

0.058243 

0.037269 

0 . 000020 

0.555717 

0.058516 

0.037370 

0 . 000030 

0.555577 

0.058782 

0.837476 

0 ,  Fin 0050 

0.555298 

0.059314 

0.037636 

0.000070 

G.5:'^5n?2 

0.059842 

0.037848 

0.000100 

0.554610 

0.060629 

0.038113 

0.000147 

0.553972 

0.061850 

0.038521 
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Table  A13.  Coordinated  failure  model,  UX  -  2.1/0.90  (Fig.  8) 


pQ  =  2.1  Mg/m^ 
=  100  kPa 

Z  =  0.1380 

e  =  0.4275 
w 


Input  variables 

S  =  0.90 
w 

=  14.7  kPa 

M 

Output 
<})q  =  0.3220 
EF  =  8.551  kJ 


K  =  5.0  GPa 


=  0.0322 
a 

EV  =  39.04  kJ 


LOADING  PATH 


p 

V 

DP/DV 

MU 

DP/DMU 

0. 

0.476190 

0. 

0. 

0, 

0.00OGO1 

0-.  476181 

-0. 105002 

0 . 000020 

0.058000 

0  j 'MNl  i  1 1'l 

O.OOOnij'i 

"0.001239 

0.015509 

0.000581 

it .  }ji  lilt  i,;ij 

0 . 4i4m'22 

"0.004553 

0.020233 

0.002093 

i'..nni;n;5n 

0.41:5.11 7 

"0 . 007663 

0.023149 

0.003496 

M  .  I'lni  !!  *:  ,11 

n.-ii,  -.m 

"0.01 1062 

0.026869 

0.005376 

!  1  .  L'lll  Ml, 'll 

0..:r.45,'n 

"0.01 7366 

0.029425 

0.007823 

i  I .  (Mil  1  liin 

o.4^4-:n4 

-0.023527 

0.0:K2271 

0.010543 

11. 4;,''ii.M: 

“0 . 03202  1 

0 . 035566 

0.014264 

1 1 .  i  iDtL  ‘i  in 

o.O'.n.inn 

"0. 140694 

0.036369 

0.065976 

n .  i.Dujni.iij 

0.4':.:UG4 

"0. 154316 

0.039301 

0.068224 

0 .  l3On(.>Li0 

0.454050 

"0.  163243 

0.042087 

0.071777 

0 . 00 1 000 

0.45.1692 

-0. 176474 

0.047232 

0.077001 

0.000000 

0.4-19053 

"0.206657 

0.058547 

0.088768 

G .  0000130 

0 . 4^i;-M36 

"0.269662 

0.076292 

0.  1 12710 

0.00,  >000 

0.4 ^G 760 

"0.552833 

8.090260 

0. 143182 

n . OD0U0O 

0.430081 

-0.443637 

0.  107212 

-  0. 176976 

0 . non 000 

0 . 4:- I'T  OiO 

"0.516684 

0.  1  10209 

0.200157 

0 .  U  1 IJOLIO 

0.4::  01 171 

"0.534266 

0.  112151 

0.206649 

0.015000 

0.420106 

"0.619965 

0.  133501 

0.234187 

0.020000 

0.413614 

"0.770209 

0.  151291 

0.281049 

O.02-4nuO 

0.4091 12 

"0.333482 

0.  163961 

0.315722 

0 . 030000 

0.403125 

-1.002115 

0.  181248 

0.347071 

0.040000 

0 . 394490 

-1 .  153080 

0.207105 

0.386755 

0. OGOonn 

0.300130 

-1.392745 

0.252705 

0.438590 

0.100000 

0.357252 

-1.743455 

0,332925 

0.498633 

0.  150000 

0. -332 194 

-1.995313 

0.433472 

0.497277 

0.200000 

0.312226 

"2.504019 

0.525148 

0.545402 

0 . 300000 

0.279110 

-3.019666 

0.706106 

0.552614 

0 . 400000 

0.262898 

-6. 163479 

0.811312 

0.950516 

Q.500G00 

0.252216 

“9.361209 

0.888028 

1.303500 

0.600000 

0.243585 

-11.586166 

0.954927 

1 . 494794 

0 . 000000 

0.230149 

-14.885425 

1.069055 

1.752425 

1 . 000000 

0.219874 

-13.464453 

1.165746 

2.068443 

1 .500000 

0.201475 

-27. 175850 

1.363522 

2.528115 

2.000000 

0.  18S567 

-33.735521 

1.525313 

3.090405 

0. 

UNLOADING  PATH 

0 . 460356 

0.033275 

0. 

0.360001 

0 . 460843 

0.033231 

0.062858 

0.0000  10 

0 . 460785 

0.033434 

0.062948 

0.000020 

0.460714 

0.033592 

0.063120 

0 . H 00030 

0 . 460644 

0.033750 

0.063301 

0.000050 

0.468503 

0.034065 

0.063572 

0.0GOO70 

0 . 460364 

0.034373 

8.063933 

0.000100 

0.460157 

0.834844 

0.864384 

0.000147 

0.459836 

0.035566 

0,065077 
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Table  A14.  Coordinated  failure  model,  SX  -  1. 7/1.0  (Fig.  8). 


Input  variables 


Iq  =  1.7  Mg/m^ 

S  =  1.0 
w 

K  =  5.0  GPa 

’  =  10  kPa 

T 

=  68.9  MPa 

M 

Output 

Z  =  0.3417 

(t)Q  =  0.5808 

(})  =  0. 

^a 

!  =  1.386 
w 

EF  =  6.612  kJ 

EV  =  33.30  kJ 

LOADING 

PATH 

P 

V 

DP/DV 

MU 

DP/DMU 

0. 

0.588235 

0. 

0. 

0. 

0.000000 

0.588234 

-0.005000 

0.000002 

0.050000 

0 . 000007 

0.588122 

-0.060403 

0.000193 

0.035527 

0.000010 

0.588071 

-0.060766 

0.000280 

0.035728 

0 .  oooo;“::0 

0.587906 

-0.060900 

0.000560 

0.035794 

0 . 000030 

0.587743 

-0.061106 

0.000338 

0.035894 

0 . 000000 

0.587417 

-0.061414 

0.001393 

0.036046 

□ .  000070' 

0.5E57094 

-0.0G1S2S 

0.001945 

0.036247 

0 . GOO  100 

0.586612 

-0.062340 

0.002767 

0.036498 

0 . 000200 

0.585042 

-0.063675 

0.005459 

0,037150 

0 . 000400 

0.582046 

-0.066753 

0.010634 

0.038646 

0 . 000600 

0.579225 

-0.070903 

0.015555 

0.040640 

0 . 00 1000 

0.574033 

-0.077113 

0.024732 

0.043588 

0 . 002600 

0.563115 

-0.09 1550 

0,044609 

0.050309 

0 . 004000 

0.546B51 

-0.122969 

0.075678 

0.064374 

0.006000 

0.534093 

-0. 167259 

0.099724 

0.083171 

0.003000 

0.521413 

-0.222552 

0.  12S155 

0. 105519 

0.003600 

0.519143 

-0.264354 

0.  133078 

0.  121879 

0.010000 

0.517700 

-0.276209 

0. 136243 

0. 126199 

0.010000 

0.502735 

-0.3341 18 

0.  170070 

0. 147831 

0.026000 

0.491549 

-0 . 44700 1 

0.  196696 

0.  137786 

G . 02 40 00 

0.404234 

-0.546760 

0.214776 

0.221242 

0.030000 

0 . 47503 1 

“0.652023 

0.238303 

0.254978 

0 . 0  "410  0  0  G 

0.462713 

-0.811816 

0.271274 

0 . 303348 

0.060060 

G . 444252 

-1.033348 

0.324102 

0.378581 

0. 100000 

0.418421 

-1.543510 

0.405846 

0.489335 

0.  100000 

0.388316 

-1.660331 

0.514335 

0.458761 

0 . 200000 

0 . 365807 

-2.221264 

0.608050 

0.536397 

0 . 300000 

0.329394 

-2.746310 

0.735810 

0.562556 

0 . 400000 

6.308441 

-4.772514 

0.907125 

0.824297 

G.5GGOG0 

0.293214 

-6.567333 

1.006164 

1.009706 

0.600000 

0.280559 

-7.902333 

1.096651 

1.195132 

0.808000 

■0.260176 

-9.811885 

1.260913 

i. 2 17569 

1 .000000 

0.243982 

-'12.350373 

1.410976 

1.332767 

1.500000 

0,213573 

-16.442434 

1.754258 

1.456531 

2 . 006000 

n.  191  106 

-22.254613 

2.078061 

1.544150 

. UNLOADING  PATH 

0 . 

0.500035 

0.000000 

0. 

0 .  no  00  OtO 

0.000034 

0 . 000003 

0 . 035643 

0.000007 

0  .  O'  *0  1 22 

0.000193 

0.035678 
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Table  A15.  Coordinated  failure  model,  SX  -  2. 1/1.0  (Fig.  7). 


Input  variables 


Pq  =  2.1  Mg/m^ 

S  =  1.0 
w 

K  =  5.0  GPa 

=  10  kPa 

=  68.9  MPa 

M 

Output 

Z  =  0,1625 

(fig  =  0.3413 

(J)  =  0. 

e  =  0.5182 
w 

EF  =  8.926  kJ 

EV  =  40.32  kJ 

LOADING  PATH 


P 
tl . 

O.0ODO00 
0 . 000007 
0.000018 
0 . 000820 
0 . 000030 
0.000050 
0 . 000070 
0. 00010^ 
0 . 000270 
0.000400 
0.000600 
0. 0010010 
0.002000 
0.004000 
0.006008 
0.009000 
0 . 009600 
0.010000 
0.015000 
0 . 020000 
G . 024008 
8 . 030008 
0.840000 
VJ  .  060000 
0.  100!  100 
0.  1501100 
0 . 200080 
0 , 300000 
0 . 4000.00 
0 . 500000 
0 .600000 
0 , 8800 80 
1.000000 
1  .500000 
2.080000 


V 

0.476190 
0.476190 
0.476134 
0.476108 
0 . 476826 
0.475945 
0.475732 
0.475621 
0 . 4753S0 
0 . 474596 
0.473895 
0.471678 
3 . 469859 
0 . 463488 
0.455886 
0 . 44853 1 
0.441076 
0.439782 
0.438949 
0.430053 
0 , 422996 
0.418156 
0.41 1781 
0.482703 
n . 38 7 049 
11 . 3i  .■  16  16 

ij.  338950 
0 . 318676 
0.285163 
0.268331 
0.257151 
0 . 243036 
0.233764 
0.222776 
0.202931 
0. 183373 


DP/DV 

0. 

-0.  105000 
-0.  121657 
“0.  121851 
-0,  122109 
-0.122503 
-0. 123092 
-0.  123879 
-0. 124862 
-0. 127408 
-0. 133273 
-0. 141 124 
-0.  152776 
-0. 179439 
-0.235798 
-0.311258 
-0.393743 
-0.463615 
-0.480250 
-0.562074 
-0.708533 
-0.826309 
-0,941272 
-1  .  101526 
-1 .346453 
-1.721694' 
-1.948095 
-2.466225 
-2.983910 
-5.958647 
-8.905133 
-10.970493 
-14.013109 
-18.202208 
-25. 195779 
-35.564763 


hU 

0. 

0.000002 

0.000119 

0.000173 

0.000345 

0.000516 

0.000858 

0.001198 

0.001704 

0.003361 

0.006543 

0.009567 

0.015203 

0.027406 

0,046558 

0.061549 

0.079611 

8.082733 

0.084843 

0. 107282 

0.125755 

0. 138787 

0. 156416 

0. 182485 

0.227772 

0.306005 

0.404898 

0.494277 

0.669888 

0.774309 

0.851791 

0.919844 

1.037060 

1.  137531 

1.346559 

1.521227 


UNLOADING  PATH 

0 .  0 . 47  6 1 90  0 . 000000 
0 . 080000  0 . 476 190  0 , 000002 
0.000087  0.476134  0.000119 


DP/DhU 

0. 

0 , 050000 

3.057925 

0,058008 

0.058117 

0.058284 

0.058535 

0.058869 

0.059286 

0.060365 

0.062840 

0.066132 

0.070382 

0.081945 

0. 104428 

0, 133413 

0.  166096 

0.  188855 

0. 194687 

0.222818 

0.278669 

0.306928 

0.340361 

0.383588 

0.441630 

0.511238 

0.585593 

0.559419 

0.569440 

0.957661 

1.230625 

1.469430 

1.706261 

1.990622 

2.392019 

2.862578 


0. 

0.057867 

0.057925 
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Table  A16.  Coordinated  failure  model,  CM  -  1.7/0.90  (Fig,  9) 


Pq  =  1.7  Mg/m^ 
=  5  MPa 


Z  =  0.2901 


e  =  1.091 
w 


P 

0 . 

n.nrinniG 

0.000030 
0 , 000050 
0 . 0!j00;oj 
0 . 000  ]  ijO 
0.000300 
n . onnp  is 

0 . 000400 
0 . OOOSOG 
0.00  1000 
0.002000 
G .  0040010 
G . 00 G 000 
0 .  Gooru-in 
0 . 00:0  04  . MO 
0.01  OiOiOiG 
Ci .  0  1500:0 
0 . 02ijni  ii:i 
0 . 02401010 
0 . 0301  iiin 

0.0401:1  in 
0 .  GGOOiOOi 
•0.  1000010 
0.  isnonn 

n .  20i;JOi.:O 
Oi .  soriioco 
0 . 400:j;:i0 
0 . 500000 
0 . 80004, 0 
G . 800000 
1 .000000 
1 . 500000* 
2.GGOOOG 


G . 0080  10 
0 . 000028 
0 . 000030 
G . 000050 
0 . 0000  ?4Zi 
0.080103 
0.000 200 
0 . GOO^To 


Input  variables 

S  =  0.90 
w 

P^  =  21.5  MPa 
M 

Output 
(t)Q  =  0.5480 
EF  =  5.939  kJ 


LOO DING  POTH 


DP/DV 

0.5:10235 

0. 

0.5001  10 

“0.085017 

0 . '  .iOOUOO 

“0.085051 

n.'5:::r383 

“0.085085 

0 . 5:37640 

“0.085136 

0.579451 

-0 . 002440 

0.571032 

“0.003563 

0 . 55-1773 

-0.006150 

0 . 553067 

“0.003796 

0 . 550633 

“0.078135 

0.54:8283 

-0.082755 

0.543835 

“0.089927 

0.534451 

“0.  106573 

0.520-426 

“0. 142600 

0 . 5  1 0060 

-0. 192940 

n .  .••r':i!:i233 

“0.255084 

n... 51 3 

“0.302120 

0.-49571-12 

“0,314674 

H. -4: 1:424 

-0.378255 

Li.47i;:26 

“0.500113 

Oi .  --  Ii  .5222 

-0.605718 

0.-15r.;:!3-l 

-0.715261 

n.-M5--i--iri 

“0.878275 

0..' .1201 124 

-1 . 147832 

0,-102541 

-1.594661 

0 . 374  1  1 3 

-1.734453 

0 . 352247 

“2.286609 

0 . 3  Ibbbo 

“2.810702 

0.296915 

-5.062453 

0 . 2::7n3F! 

-7. 103901 

0.271202 

“3.593845 

n .  252577 

-10,733134 

0.2378S1 

-13.609175 

0.2  1051 1 

-13.268455 

0.  190463 

-24.940076 

UNLOADING  PATH 

8.555933 
0.555858 
0.555717 
0,555577 
01 . 55529o 
0.555022 
0.554610 
3.553265 
0.553067 


K  =  5.0  GPa 


MU 

(f)  =  0.0548 

a 

EV  =  31.77  kJ 

DP/DMU 

0. 

■  0. 

0.000200 

0 . 050000 

0.000400 

0.050000 

0.000600 

0.050000 

0.001000 

0.050000 

0,015159 

0.001413 

0.030126 

0.002004 

0.060318 

0.003312 

0.063587 

0.004538 

0.068160 

0.040456 

0.072369 

0.042473 

0.081644 

0.045584 

0.  100634 

0.052659 

0. 130295 

0.067427 

0.  153266 

0.087067 

0.  180486 

0.  110216 

0.  185209 

0.  127021 

0.  188253 

0. 131434 

0.220852 

0.  153378 

0.246721 

0.  193280 

0.264418 

0.226023 

0.287636 

8.258424 

0.320543 

0.303333 

0.374303 

0.372059 

0.459856 

0.467549 

0.572346 

0,444483 

0.669953 

0.512260 

0,857575 

0.532985 

0.981156 

0.809183 

1.079758 

1.014183 

1 . 168992 

1.  120644 

1.328936 

1,250444 

1.472814 

1.390061 

1.794316 

1.555199 

2.088445 

1.699939 

0.057980 

0. 

0.058243 

0.037264 

0 . 0585 16 

0.037370 

0.058782 

0.037476 

•0.059314 

0.037636 

0.059842 

0.037848 

0.060629 

0,033113 

0.063207 

0.038799 

0.063587 

0.039410 
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Table  A17.  Coordinated  failure  model,  CM  -  2.1/0.90  (Fig.  9) 


Input  variables 


pQ  =  2.1  Mg/m^ 
=  5  MPa 


Z  =  0.1380 


e  =  0.4275 


P 

0. 

0 . 0000  10 
0 . 000020 
0 . 000020 
0 . 000050 
0.000070 
0.000100 
0.080200 
0.000215 
0 . 000400 
0 . 0OOGO0 
0 . 001000 
0 . 002000 
0 . 004000 
0 . 005000 
0.005000 
0.005500 
0 . 0  1 0000 
0 . 015000 
0 . 020000 
0 . 0240UO 
0 . 03GOOO 

4.!i4;:5l;0 

ii.  inn: 00 

0.  lo: non 

i  '  .  i  :-  '’Oj 
0.:vM:i7:n 

0 .  nooLiOO 
\) .  5i  in  ion 
0 . 1*0!  If  sni 

n  .on.noo 
i  .nnnnoo 
1  .  OilOliOO 
2.000000’ 


n . 

I  ].  00  on  10  tj 
n.n:ion20  O 
0.000020  G 
0.0!  no  11 50  0 
0.000070  G 
0.000100  0 
0.000200  0 
0.000215  0 


S  =  0.90 
w 

P  -  21.5  MPa 
M 


Output 
(j)^  =  0.3220 
EF  =  8.552  kJ 


LGhDING  path 


y 

DP/DV 

0.4?bl90 

0. 

0.475095 

-0. 105021 

0,476000 

-0. 105063 

0 . 475905 

-0. 105105 

G. 4757 15 

-0. 105168 

0.472059 

-0.005471 

U .  4'3UOb7 

-0.007515 

0.4G0214 

-0.012734 

0.459380 

-0.017971 

0.^53184 

-0. 154653 

0.456958 

-0. 163248 

0.454692 

-0. 176474 

0 . 449853 

-0.206657 

0.442436 

-0.269662 

0 . 436768 

-0.352833 

0.430031 

-0.448637 

0 . 428920 

-0.516634 

0.428171 

-0.534266 

0.420106 

-0.619965 

0.413614 

-0.770209 

0.409112 

-0.883482 

0 . 403  125 

-1.002115 

0.  O  14490 

-1  .  158088 

U  .  3:.:0  1  30 

-1.392745 

0.357252 

-1.743455 

n. 332 194 

-1 .995313 

0  .  :•  1  22? 

-2. 5040 19 

0.2391  in 

-3 . 0  19666 

0 . 252898 

-6 . 1 63479 

1^252216 

-3.361209 

0.24 3585 

-1  1.586166 

0.230149 

-14.885425 

0.2I98?4 

-19.464453 

0.2!:U475 

-27. 175850 

0-  188567 

-38.735521 

UNLOAD  I HG  PATH 

0 . 4 50 05 G 
450705 
450714 
4GG644 
4G05G3 
450354 
4G0  157 
459400 
459330 


K  =  5.0  GPa 


d)  -  0.0322 

^a 

EV  =  39.04  kJ 


MU 

DP/DMU 

0. 

0. 

0.000200 

0.050000 

0 . 000400 

0.050000 

0.000600 

0.050000 

0.001000 

0.050000 

0.008751 

0.0025S0 

0.017355 

0.003487 

0.034714 

0.005760 

0 , 036594 

0.007979 

0.039301 

0.068360 

0.042037 

0.071777 

0.047232 

0.077001 

0.058547 

0.088768 

0.076292 

0. 112710 

0.090260 

0. 143182 

0.  107212 

0.  176976 

0. 1 10209 

0 . 200157 

0.  1  12151 

0.206049 

0.  133501 

0.234187 

0.  151291 

0.281049 

0. 163961 

0.315722 

0.  131248 

0.347071 

0.207105 

0.336755 

0.252705 

0.438590 

0.332925 

0.493633 

0.433472 

0.497277 

0.525148 

0.545402 

0 . 706 106 

0.552614 

0.311312 

0.950516 

0.833028 

1.303500 

0.954927 

1.494794 

1.069055 

1.752425 

1. 165746 

2.068443 

1.363522 

2.523115 

1.525313 

3.090405 

.033275 

0. 

,...033434 

0.062939 

0 . 033592 

0.063120 

0.033750 

0.063301 

0.034065 

0.063572 

0.034378 

0 , 063933 

0.034844 

0.064334 

0 . 036369 

0.065550 

0.036594 

0.066587 
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Table  A18.  Coordinated  failure  model,  CX  -  1.7/0.90  (Fig.  9) 


Z  = 

e  = 
w 


Input  variables 


!•?  Mg/m^ 

S  =  0.90 
w 

K  = 

=  5.0  GPa 

5  MPa 

=  25.1  MPa 

M 

Output 

0.2901 

())q  =  0.5480 

K 

=  0.0548 

1.091 

EF  =  5.939  kJ 

EV 

=  31.77  kJ 

LOODING  PATH 

P 

y 

DP/PV 

MU 

DP/DMU 

0. 

0 . 588235 

0. 

0. 

0. 

6.000010 

0.5881  is 

-0.085017 

0.000200 

0 . 050000 

0 . 000020 

0 . 5:0iJ0O8 

“0 , 0S505 1 

0.000400 

0.050000 

0.000030 

0, 51. i  788  3 

-0.085035 

0.000600 

0 . 050000 

0 .  uoni  jSO 

0 . 587648 

-0 . 085 136 

0.001000 

0 . 050000 

0 . 000070 

0.580186 

-0.002630 

0.013873 

0.001554 

0 . 000100 

0.572505 

-0.003906 

0.027476 

0.002205 

0 . 0G02Otn 

0 . 55f7567 

-0.006694 

0 . 055004 

0.003633 

0.000251 

0.552596 

-0.010260 

0.064495 

0.005374 

0 .  ijooono 

0 . 550693 

-0.078574 

8.068166 

0.040649 

0 .  Of!  Ob  on 

0 . 548283 

-0 . 082755 

0.072869 

0 . 042478 

0 .  or;  1  ni'in 

0 . 543035 

-0.089927 

0.081644 

8.045584 

f! ,  Of  i.  -i  If  !!;i 

1 1 . 5  7^15  I 

-n. 106573 

0.  100634 

0.052659 

0 . 1  ii  i.  ii  tr.i.i 

0  .  0:'G 

-0. 142600 

0. 130295 

0,067427 

l.l  .  1  iOi ;!  IIJ 

0.5  lOi.ir.iO 

-0. 192340 

0.  153266 

0.037067 

0  .  i'O’.OOt  10 

n .  4'Xi239 

-0.255084 

0. 130486 

0.  1102‘16 

1  1  .  i  11  i-  !|  ,1  il  l 

0.456313 

-0.302120 

0. 185209 

0. 127021 

11.0  101  100 

0 . 4’'i5n42 

-0.314674 

0.  1CS253 

0.  131434 

U  .  0  15000 

0 . 4:;::  1024 

-0.378255 

0.220852 

0. 153378 

0 . 020000 

0.471026 

-0.5001 13 

0.246721 

0. 193288 

0 . 024000 

0 . 465222 

-0.605713 

0.264418 

0.226028 

0 . 030000 

0 . 456034 

-0.715261 

0.287636 

0.253424 

0 .  n.  ionnri 

0 . 445448 

-0.878275 

0.320548^ 

0.303833 

n .  nr.  0000 

n .  4; 33024 

-1. 147832 

0.374303 

8.372059 

0.  iOOi  ii  in 

n.4!.i254i 

-1.594661 

0.459856 

0.467549 

0.  15001:0 

0.374113 

-1 .734453 

0.572346 

0 . 444433 

0.200000 

0.352247 

-2.286609 

0.669953 

0.512260 

n.sonocin 

0  .  -3  J.  b  6  b  cl 

-2.310702 

0.857575 

0.532985 

0 . 400000 

0.296915 

-5.062453 

0.931156 

0.809183 

0.500000 

n . 202033 

-7.  103901 

1.079758 

1.014183 

0 . 600000 

0.271202 

-8 , 593345 

1. 168932 

1. 120644 

0.800000 

0.252577 

-10.738134 

1.328936 

1.250444 

1 . 000000 

0 . 23700 1 

-13.609175 

1.472814 

1.330061 

1 . 500000 

0.210511 

-13.268455 

1.794316 

1.555199 

2 . 0801100 

0.  190463 

-24.940076 

2.088445 

1.693939 

UNLOADING  PATH 

0. 

0.555999 

0.057930 

0. 

0 . 0000  10 

0 . Obbobo 

0.058248 

0.037264 

0 . 0060O:j 

0.555717 

0.058516 

0.037370 

0 . 000000 

0.555577 

0.058732 

0.037476 

0.000050 

8.555293 

0.059314 

0.037636 

0.000070 

8.000100 

-4  qcTirr^O-p 

0 . 5 5 46 7 8 

0.059842 

0.060629 

0.037348 

0.038113 

0.000200 

0.553265 

0 . 063207 

0.038799 

0.000251 

0.552596 

0.064495 

0.039599 
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Table  A19.  Coordinated  failure  model,  CX  -  2.1/0.90  (Fig.  9) 


Pq  =  2.1  Mg/m^ 
=  5  MPa 

Z  =  0.1380 
e  =  0.4275 


Input  variables 

S  =  0.90 
w 

-  25.1  MPa 

M 

Output 
(f)Q  =  0.3220 
EF  =  8.552  kJ 


K  =  5.0  GPa 


(f)  =  0.0322 

a 

EV  =  39.04  kJ 


LOADING  FRTH 


P 

0. 

0.000010 
0 . 800020 
0.000030 
0.000050 
0.000070 
0.000100 
0 . 000200 
0.000251 
0.000400 
0.000606 
0.001000 
0 . 002000 
0.004000 
0 .  EiObOrjO 
6.00SDO0 
0.000600 
0,010000 
0.015000 
0.020000 
0.024000 
0.030000 
0 . 040000 
0.060000 
0;  100000 
0. 156000 
0 . 200000 
0.300860 

0.400000 

0.500000 

O.bnnnnn 

.  WOO 


V 

0.476190 
0.476095 
0 . 476000 
0.475305 
0.475715 
0 . 472405 
0 . 4b 6  767 
0.461572 
0.459142 
0.453134 
0 . 456958 
0.454692 
0.449853 
0.442436 
0.436768 
0 . 43003 1 
0 . 428920 
0.428171 
6.420106 
0.413614 
0.409112 
0.493125 
0.394490 
0.330130 
0.357252 
0.332194 

0 . 4 1 2226 

0.279110 

0.262898 

0.252216 

0.245535 


V  .  i  517567 


DP/DV 

0. 

-0.105021 
-0. 105063 
-0. 105105 
-0.105163 
-0.006042 
-0.008243 
-0.013398 
-0.020986 
-0. 155476 
-0.163243 
-0. 176474 
-0.206657 
-0.269662 
—0 . 352333 
-0.443637 
-0.516634 
-0 . 534266 
-0.619965 
-0.770209 
-0.88S482 
-1.002115 
-1. 158088 
-1.392745 
-1.743455 
-1.995318 
-2.504019 
-3.019666 
-6. 168479 
•'-9.361209 
-1 1.586166 
-14.335425 
-19. 464453 
-27. 175350 
-38.735521 


MU 

0. 

0.000200 
0 . 000400 
0 .000600 
0.001000 
0.008014 
0.015835 
0.031671 
0.037131 
0.039301 
0.042087 
0.047282 
0.058547 
0.076292 
0.030260 
0. 107212 
0. 1 10209 
0.  112151 
0.  133501 
0. 151291 
0. 163961 
0.181248 
0.207105 
0.252705 
0.332925 
0.433472 
0.525148 
0.706106 
0.811312 
0.888028 
0.354927 
'1.069055 
1 . 165746 
1.363522 
1 .525313 


DP/DMU 

0. 

0 . 050000 

0.050000 

0.050000 

0.050000 

0.002851 

0.003836 

0.006315 

0 , 009340 

0.068686 

0.071777 

0.077001 

0.088768 

0.112710 

0. 143182 

0. 176976 

0.200157 

0.206049 

0.234187 

0.281049 

0.315722 

8.347071 

0.386755 

0.438590 

0.498633 

0.497277 

0.545402 

8.552614 

0.950516 

1.303500 

1 . 494794 

1.752425 

2.868443 

2.528115 

3.090405 


i:ill  UmDING  path 


1 1 . 

0.033275 

0. 

1.; . ;  .1  iiU  1 1 0 

11 . 

0.033434 

0.062939 

i  ; ,  ni  !i  H  t.M:i 

0. 

.;i.riri4 

0.033592 

0.063120 

ft .  I,:!  -i 

0. 

•  !’..i;n;4.) 

0.033750 

0.063301 

.  i  .M 

Ti . 

U . 034065 

0.063572 

I  t .  Onij:  ii’0 

1 1 . 

il-lKMLl 

0.034378 

0.063933 

I'i .  iJUt  i  11  111 

0  . 

0.034844 

0.064384 

0 .  Ui'K i'nitj 

0  . 

.i:.,9'180 

0.036369 

0.065550 

0.1100?':^ 

0  . 

459142 

0.037131 

0.066908 
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Table  A20.  Coordinated  failure  model,  RM  -  1,7/0.90  (Fig.  10) 


Input  variables 


)q  =  1.7  Mg/m^ 

S  =  0,90 
w 

K  = 

=  5.0  GPa 

=  8  MPa 

=  34.5  MPa 

M 

Output 

Z  =  0.2901 

<()p  =  0.5480 

=  0.0548 

!  =  1.0913 
w 

EF  =  5.940  kJ 

EV 

=  31.77  kJ 

LOOniMG  PATH 

P 

V 

DP/DV 

MU 

DP/DMU 

0. 

0.30B235 

0. 

0. 

0. 

o.oofitun 

n.ORRUO 

-0.085017 

0.000200 

0.050000 

U  .01.11  li  t i-'iJ 

0 .  l.-:.,!0nLiR 

-0,085051 

0 . 000400 

0.050000 

0. 000;  130 

i:\  r .  1  1  *vri  0-7 

-8.085085 

0.000600 

0.050000 

G.OOUOSO 

0 . 3 0/6 48 

-G.G35136 

0.001000 

0.050000 

0 .  onriopR 

0.30?^:n5 

-0.085204 

0.001400 

0.050000 

D.OOO.iOR 

R . 3 87 29b 

-0.085255 

0.001600 

0.050000 

O.nor.jOO 

o.ooiai? 

-0.003522 

0:011379 

0.002045 

o.ooopon 

0 . 564812 

-0.005951 

0.041470 

0 . 003323 

Cl  . 

Ij4j51389 

-0.010802 

0.066824 

0.005719 

U  ,  LiOLMUIJ 

0 . 550699 

-0.079712 

0.068160 

0.041 148 

0 . nnoono 

0.548283 

-0.082755 

0.072869 

0.042478 

0 . GG 1 nnn 

0.543835 

-0.089927 

0.081644 

0 . 045584 

0 .  O02OJU0 

0 . 53445 1 

-0.  106573 

0.  100634 

0.052659 

0 . 00«^]nn0 

0.520426 

-0. 142600 

0.130295 

0.067427 

0 .  nnonriR 

0.51 0060 

-0.  192940 

0. 153266 

0.087067 

0 .  nno’ion 

0.498299 

-0.255634 

0 . 1 80486 

8. 1  10216 

o.Gii'j.  tin 

0.496313 

-0.302120 

0.  185209 

0 . 12702 1 

c .  n  iG.ujR 

0 . 495042 

-0.314674 

0.  188253 

0. 131434 

0.01G0,..J 

0.481824 

-0.378255 

0.220852 

0. 153378 

0 . 020000 

0.47 1826 

-0.500113 

0.246721 

0. 193280 

G . 024000 

0 . 465222 

-0.605718 

0.264418 

0.226028 

G . 030000 

0 . 456834 

-0.715261 

0.287636 

0.258424 

0.040000 

0.445448 

-0.878275 

0.320548 

0.303833 

0.060000 

0.428024 

“1. 147882 

0.374303 

0.372059 

0. 100000 

0 . 402941 

-1.594661 

0.459856 

0.467549 

0. 150000 

-  0.3741 13 

-1.734453 

0.572346 

0 . 444433 

0 . 20000(0 

P  ■■■;,  tl, 

-2.236609 

0.669953 

0.512260 

0.3000(00 

0  .  6  1  bbSo 

-2.810702 

0.857575 

0.532985 

0.4000010 

0.296915 

-5.062453 

0.981156 

0.809183 

0.500000 

0.282838 

-7.  103901 

1.079753 

1.014183 

G .  600000 

0 . 271202 

-8.593845 

1 .  168992 

1 .  120644 

D. oOGOOO 

0 . 2525177 

-10.738134 

1.328936 

1.250444 

1 . 000300 

0.237881 

-13.609175 

1.472814 

1.330061 

1.500000 

0.210511 

-18.268455 

1.794316 

1.555199 

2.000000' 

0.  190463 

-24.940076 

2.088445 

1.693939 

UHLuin-'ING  POTH 

0 . 05 7 930 

0. 

O/Z/tOOR 

0.058248 

0.037264 

0.000717 

0.058516 

0 . 037370 

6. 00' -.0  7  7 

0.0587S2 

0.037476 

0 . 0001.100 

0.00Z.290 

0.059314 

0.037636 

i' .  000070 

0.0':!0D22 

0.059842 

0.037848 

0.(700000 

0.004384 

0.060105 

0.038037 

0  .  L'iJu  j  (00 

0.004010 

0 . 060629 

0.038166 

0 . UR 02 00 

0  0  Z.  ‘s'>2b0i 

0.063207 

0.038799 

0 . 000345 

0.001089 

0.066824 

0.040091 
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Table  A21.  Coordinated  failure  model,  RM  -  2.1/0.90  (Fig.  10) 


=  2.1  Mg/m^ 
=  8  MPa 

Z  =  0.1380 
e  =  0.4275 


Input  variables 

S  =  0.90 
w 

Pj^  =  34.5  MPa 

Output 
(j)^  =  0.3220 
EF  =  8.552  kJ 


K  =  5.0  GPa 


(()  =  0.0322 

EV  =  39.04  kJ 


LOODING  PATH 

p 

DP/DV 

MU 

DP/DMU 

0. 

0 . 47  6 1 50 

0. 

0. 

0. 

n. 0000 10 

0.476095 

-0. 105021 

0.000200 

0.050000 

0.000020 

0.476000 

-0. 105063 

0.008400 

0.050000 

0.000030 

0 . 475305 

-0. 105105 

0.000600 

0.050000 

0 . 000000 

Pj  .  4757  1 5 

-0. 105163 

0 . 00 1000 

0 . 050000 

0.000070 

0.475525 

-0. 105252 

0.061400 

0.050000 

n .  orifinon 

0 . 475430 

-0. 105315 

0.001600 

0 . 050000 

11  .VII  10  ]  00 

0.473065 

-0.003453 

0.006606 

0.003995 

o.RoiiOon 

0. -165067 

-0.F1 12502 

0.023919 

0.005776 

0 . 000340 

0 . 450533 

-0.022191 

0.033510 

0.009938 

0 . 000400 

0.453104 

-0. 157593 

0.039301 

0.069531 

\i ,  onoron 

0 . 456050 

-0. 163243 

0.042087 

0.071777 

1 1 . 1 11  M  000 

0.454692 

-0. 176474 

0.047232 

0.077001 

n . no? 000 

0.443B53 

-0.206657 

0.058547 

0.088768 

0.004000 

0 . 44243G 

-0.2696G2 

0.076292 

0. 112710 

0 .  onr.rii'!::i 

0 . 43G7C8 

-0.352833 

0.090260 

0. 143182 

0  .  OGOOriei 

0 . 43008 1 

-0.443637 

0. 107212 

0.  176976 

0 . 000600 

0.428920 

-0.516684 

0. 110209 

0.200157 

n :  0  iCionn 

0.428171 

-0.534266 

‘0.112151 

0.206049 

0.015000 

0 . 420 1 06 

-0.619965 

0. 133501 

0.234187 

n. 020000 

0.433614 

-0.770209 

0.151291 

0,281049 

0.024000 

0.409112 

-0.838432 

0.  163961 

0.315722 

0 . 030000 

0.403125 

-1.002115 

0. 181248 

0.347071 

0 . 0411000 

0 . 394490 

-1  .  153033 

0.207105 

0.3S6755 

0 . 000000 

0 . 3  8  Li  1 30 

“1 .392745 

0 . 252705 

0.438590 

D. in no 00 

Li  .  3 57 2 5? 

-1.743455 

0.332925 

0.493633 

0  .  150 LI 00 

l;!332'i94 

-1.995318 

0.433472 

0.497277 

G . 200000 

0 . 3  1 222)3 

-2.504019 

0.525148 

0.545402 

0 . 300000 

0!279T'i0 

-3.019666 

0.706106 

0.552614 

O.40Gnn0 

0.262893 

“6 . 16S4i-9 

0.81 1312 

0.950516 

G . 500000 

0.252216 

-9.361209 

0,833028 

1.363500 

D.buGGGG 

0,243535 

-1  1.536166 

0.954927 

1 .494794 

0.SCGC-GG 

0.3:30149 

-14.835425 

1.069055 

1.752425 

1  .  1-1000;  iO 

0.219874 

“19.464453 

1, 165746 

2.068443 

1.500000 

0.201475 

-27. 175350 

1.363522 

2.528115 

2 . OGG000 

0. 18B567 

-33.735521 

1.525313 

3.090405 

6 . 

UNLOADING  PATH 

0 .  *^60856 

0.033275 

0. 

0 . 0000  10 

0 . 460785 

0.033434 

0.062939 

0 . 000020 

0.460714 

0.033592 

0.063120 

0 .000060 

0 . 400644 

0.033750 

0.063301 

0.000059 

0 . 460503 

0.034065 

0.063572 

0.000070 

:G. 460364 

0.034378 

0.063933 

0 . 060080 

.0 . 468295 

6.034533 

0.064204 

0.000100 

G . 460 15? 

6.034844 

0.064474 

0 . 000200 

0.45S480 

0.036369 

0.065550 

1..!:;  :i  :  L  . 

0.038510 

0.067741 

-67- 


Table  A22.  Coordinated  failure  model,  RX  -  1.7/0.90  (Fig.  10) 


Z  = 

e  = 
w 


Input  variables 


1.7  Mg/m^ 

S  =0,90 
w 

K  = 

=  5.0  GPa 

8  MPa 

=  60  MPa 

M 

Output 

0.2901 

(|)q  =  0.5480 

=  0.0548 

1.091 

EF  =  5.941  kJ 

EV 

=  31.78  1 

LOADING  PATH 


P 

V 

DP/DV 

MU 

DP/DMU 

0  . 

0.588235 

6. 

0, 

0. 

0.000010 

0 . 5  0  8  1  IS 

-0.085017 

e. 000200 

0.050000 

0.000020 

0.588000 

-0.085051 

0 , 000400 

0.050000 

!; . \h-\ 

0 . ' i.'pp:-; 

-n . 085085 

0.000600 

0.050000 

i  ;  .  i  1  :!  u  )'  ,i  i 

m.p::2i,44 

-  0 . 085 1 36 

0.001000 

0.050000 

PI  M 

P.‘':;r4  13 

-0.085204 

0.001400 

0.050000 

i . .  1  ;i  )i  M 

!4  P:  ..'Pbp 

-0.005255 

0% 00 1600 

0.050000 

IPO 

0  .'.i:<2iL.'2 

~0 . O046S0 

0.003942 

0.002724 

! : .  p;  ;r;  -m  in 

0 .  ‘../nxi  16 

-0 . 007864 

0.031438 

0.004445 

n.  PPL-’!  in 

i..i .  55!-.;.M0 

-0.014852 

0.056381 

0.008018 

0  .  DP  OP!  IP 

O.P  i0.4-!3 

-0.023371 

0.072869 

0.012130 

p .  Pl.i  1  pno 

0.543035 

-0.08992? 

0.081644 

0.045584 

G .  pru'iji  iO 

O.LPyMOl 

-0.  106573 

0.  100634 

0.052659 

tJ .  P0-ini':O 

0.520426 

-0. 142600 

0.  130295 

0.067427 

0 . 000000 

0.510060 

-0. 192940 

0.  153266 

0.087067 

n . 000000 

0.450209 

-0.255084 

0. 1304S6 

0.  1 10216 

0 . GOObUn 

0.436313 

-0.302120 

0.  185209 

0.  127021 

0 . 010000 

,0.435042 

-0.314674 

0.  188253 

0.  131434 

0.015000 

10.401824 

-0.378255 

0.220852 

0. 153378 

0.020000 

41.^171826 

-0.500113 

0.246721 

0.  193280 

0 . 024000 

0.465222 

-0.605718 

0.264418 

0.226028 

D  .  0:.  01-1 00 

0 .  ^15Go.34 

“0.715261 

0.287636 

0.258424 

0 . 040000 

0 . 445448 

-0.878275 

0.320548 

0.303833 

0 . 060000 

0.428024 

-1.147882 

0.374303 

0.372059 

0.  100000 

0.402941 

-1.594661 

0.459856 

0.467549 

0.  150000 

0.3741 13 

“1.734453 

0.572346 

0.444433 

G . 200000 

0 . 352247 

-2.286609 

0.669953 

0.512260 

0 . 300000 

0.316668 

-2.810702 

0.857575 

0.532985 

0 . 400000 

0.296915 

-5.062453 

0.931156 

0.803183 

0.5 00 non 

0 . 282838 

-7. 103901 

1.079758 

1.014183 

0 . 600000 

0.271202 

-8.593845 

1.  168992 

1. 120644 

0 .  oGOijiJO 

0 . 2525?'*? 

“10.738134 

1.328336 

1.250444 

1 .  ODOCii'JO 

G  .  23  7  b' 3  1 

-13.609175 

1.472814 

1.390061 

1 .500000 

0.21051  1 
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Table  A23.\  Coordinated  failure  model,  RX  -  2.1/0.90  (Fig.  10). 
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Table  A24.  Coordinated  failure  model,  CX  -  1.7/0.90  (2.5) (Fig.  9) 
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Table  A25*  Coordinated  failure  model,  CX  -  1.7/0.90  (7.5) (Fig.  9) 
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